A© 8 40 51 9 


AFWL-TR-68-46 


AFWL-TR- 

68-46 


EXPONENTIALLY  DECAYING  PRESSURE  PULSE 


MOVING  WITH  SUPERSEISMIC  VELOCITY 
ON  THE  SURFACE  OF  A  HALF  SPACE  OF 
VON  MISES  ELASTO-PLASTIC  MATERIAL 

Hans  H.  Bleich  Alva  Matthews 

Paul  Weidlinger,  Consulting  Engineer  P 
New  York,  New  York  10017 

Contract  F29601-67-C-0091 

TECHNICAL  REPORT  NO.  AFWL-TR-68-46 

August  1968 

AIR  FORCE  WEAPONS  LABORATORY 

Research  and  Technology  Division 
Air  Force  Systems  Command 
Kirtland  Air  Force  Base 
New  Mexico 


This  document  Is  subject  to  special  export  controls  and  each  transmittal 
to  foreign  goverpmeni®, ^  or  foreign  nationals  may  be  made  only  with  prior 
approval  of  AFKL  (WLDCJ  ,  Kirtland  AFB,  KM,  87117. 


AFWL-TR-68-46 


EXPONENTIALLY  DECAYING  PRESSURE  PULSE 
MOVING  WITH  SUPERSEISMIC  VELOCITY 
ON  THE  SURFACE  OF  A  HALF  SPACE 
OF  VON  MISES  ELASTO-PLASTIC  MATERIAL 

Hans  H.  Bleich  Alva  Matthews 

Paul  Weidlinger ,  Consulting  Engineer 
New  York,  New  York  10017 
Contract  F29601-67-C-0091 


TECHNICAL  REPORT  NO.  AFWL-TR-68-46 


This  document  is  subject  to  special  export  controls 
and  each  transmittal  to  foreign  governments  or 
foreign  nationals  may  be  made  only  with  prior 
approval  of  AFWL  (WLDC) ,  Kirtland  AFB,  NMex  87117. 
Distribution  is  limited  because  of  the  technology 
discussed  in  the  report. 


AFWL-TR-68-46 


FOREWORD 


This  report  was  prepared  by  Paul  Weidlinger,  Consulting  Engineer,  New  York, 
N.  Y.  under  Contract  F29601-67-C-0091*  The  research  was  performed  under 
Program  Element  6.16.46.01.41,  Project  5710,  Subtask  RSS2144,  and  was  funded 
by  the  Defense  Atomic  Support  Agency  (DASA) . 

Inclusive  dates  of  research  were  October  1967  to  May  1968.  The  report  was 
submitted  18  June  1968  by  the  AFWL  Project  Officer,  Dr.  Henry  F.  Cooper,  Jr. 
(WLDC) . 

This  report  has  been  reviewed  and  is  approved. 


DR.  HENRY  F.  COOPER,  JR. 
Project  Officer 


ROBERT  E.  CRAWFORZK 
Lt  Col  USAV 

Chief,  Civil  Engineering  Branch 


^  GEORGE  C.  DARBY  JR/  ' 
Colonel  USAF 

Chief ,  Development  Division 


ii 


i 


ABSTRACT 


(Distribution  Limitation  Statement  2) 


An  approximate  solution  is  given  for  the  effect  of  an 
exponentially  decaying  pressure  pulse  traveling  with  super- 
seismic  velocity  on  the  surface  of  a  half-space.  The  half¬ 
space  is  an  e las t i c- p las t ic  material  of  the  von  Mises  type. 
The  effect  of  a  step  wave  for  this  geometry  and  medium  was 
treated  previously.  For  that  case,  the  peak  pressures  do 
not  decrease  with  increase  in  depth,  while  such  a  decrease 
is  obtained  for  a  decaying  surface  load.  The  prime  purpose 
of  this  investigation  is  to  determine  the  magnitude  of  this 
attenuation.  The  approximate  solutions  obtained  are  valid 
for  a  limited  distance  behind  the  wave  front,  and  are  tabu¬ 
lated  for  different  sets  of  parameters  pertaining  to  the 
material  and  velocity.  The  tabulated  results  show  that  the 
peak  pressures  in  the  case  of  the  decaying  surface  load  do 
decrease  with  depth,  but  that  the  decrease  is  less  than  one 
might  intuitively  expect.  On  the  other  hand  the  attenuation 
is  in  general  larger  than  that  encountered  in  the  similar 
problem  of  an  elasto-plastic  material  of  the  Coulomb  type. 
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SECTION  I 


INTRODUCTION 


In  a  preceding  analysis,  Refs.  [1]  and  [2]  ,  the  effect 

of  a  step  pressure,  Fig.  1,  progressing  with  a  superseismic 
velocity  V  on  the  surface  of  a  half-space  has  been  studied 
for  an  elastic-plastic  material  subject  to  the  von  Mises  yield 
condition 


J2  -  ic  =  0 


(1) 


where  k  is  the  yield  stress  in  shear  and  is  the  invariant 


J  2  =  2  SijSij 


(2) 


The  present  report  considers  the  more  general  problem 
of  a  decaying  pressure  pul  p(x  -  Vt),  Fig.  2,  moving 
with  a  superseismic,  constant  velocity  V.  The  problem  under 
study  here  is  more  realietic  than  the  one  treated  in  Refs.  [1] 
and  [2],  but  it  is  also  considerably  more  complex.  Because 
of  the  complexity  only  an  approximate  solution  will  be  derived, 
which  is  valid  for  a  limited  distance  from  the  moving  front. 


[1]  Bleich,  H.H.  and  Matthews,  A.T.,  "Step  Load  Moving 
with  Superseismic  Velocity  on  the  Surface  of  an 
Elastic-Plastic  Half-Space",  Int.  J.  Solids  and 
Structures,  3,  819-852,  1967.  (Also  Office  of  Naval 
Research,  Tech.  Rpt.  38,  Contract  Nonr- 266 (86 ) , 

December  1965 . ) 

[2]  Matthews,  A.T.  and  Bleich,  H.H.,  "Stresses  in  an 
Elastic-Plastic  Half-Space  Due  to  a  Superseismic 
Step  Load",  TBL,  Ballistics  Research  Laboratory, 

Aberdeen,  Md . ,  Tech.  Rpt.  4,  Contract  DA-30-069-AMC-8 (R ) , 
March  1966. 
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FIG.  i 


Step  Pressure  moving  over  Half-Space 
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FIG.  2 

Decaying  Pressure  moving  over  Half-Space 


2 


The  approximation  emoloyed  is  somewhat  similar  to  the  one 

* 

utilized  and  discussed  in  Ref.  [3]  for  a  material  subject 
to  a  different  yield  condition. 

Because  the  velocity  V  of  the  load  is  constant,  the 
resulting  stress  and  velocity  fields  remain  unchanged  if  con¬ 
sidered  in  a  coordinate  system  moving  with  velocity  V.  In 
other  words,  Refs.  [1],  [2]  and  [3]  and  the  present  paper 
consider  steady-state  solutions,  which  may  be  used  as  approxi¬ 
mations  to  actual  situations  when  the  velocity  V  varies 
gradually . 

It  was  found  in  Refs.  [1]  and  [2]  that  the  response  to 
the  step  pressure  changes  in  character  depending  on  the  range 
in  which  the  values  of  Poisson  ratio  v,  of  the  velocity  ratio 
V/Cp  ,  and  of  the  intensity  of  the  surface  step  load  p^  are 
situated.  As  introduction  to  the  problem  with  a  decay  of 
pressure,  some  aspects  of  the  results  of  the  analysis  for  the 
step  pressure  are  restated  in  the  following. 

Defining  the  location  of  a  field  point  by  the  polar 
coordinates  <j>  and  r,  Fig.  1,  where  the  origin  of  the  coordi¬ 
nate  system  lies  at  and  moves  with  the  front  of  the  surface 
load,  it  was  found  that  the  stresses  and  velocities  ate 
functions  of  the  angle  <{>  only.  For  a  step  pressure  there  is 


[3]  Bleich,  H.H.  and  Matthews,  A.T.,  "Exponentially 
Decaying  Pressure  Pulse  Moving  with  Supe^-seismic 
Velocity  on  the  Surface  of  a  Half-Space  of  Granular 
Material",  Tech.  Rpt.  AFWL-TR-67-21 ,  Contract 
AF  29  (  601 1 -  7  08  2 ,  July  1967. 
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therefore  no  dependence  of  stresses  or  velocities  on  the 
radius  r,  so  that  there  is  no  attenuation  with  depth.  It  was 
found,  further,  that  for  sufficiently  large  values  of  the 
surface  load,  pQ  >  p^  ,  the  solution  has  the  character  of 
Fig.  3,  where  the  value  of  p^  depends  on  the  parameters  v 
and  V/Cp  . 

The  configuration  shown  in  Fig.  3  shows  two  discon¬ 
tinuities,  a  P-front  and  an  S-front,  moving  with  the  veloci¬ 
ties  of  elastic  P-  and  S-waves 


2  =  2  ( 1  -  V )  G 
CP  1  -  2V  p 


(3) 


respectively.  The  locations  of  the  fronts  are  defined  by 
the  related  angles 


The  changes  in  stress  and  velocity  at  these  fronts  are 
entirely  elastic.  There  are  four  "neutral"  regions  of 
uniform  stress  (and  uniform  velocity),  i.e.,  regions  where 
the  yield  condition,  Eq.  (1),  is  just  satisfied,  but  no  yield 
occurs  at  the  particular  instant.  There  are  two  plastic 
regions,  <{>^  to  4> ^  >  and  from  ^  to  <j>^  ,  where  plastic  defor¬ 
mations  occur  at  the  particular  instant, and  where  the  stresses 
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*  NEUTRAL  REGIONS  OF 
UNIFORM  STRESS 
SATISFYING  EQ.l 


CONFIGURATION  FOR  STEP  PRESSURES  P0^  PL 
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(and  velocities)  vary  as  functions  of  the  location  <{> .  One 
of  the  plastic  regions  is  always  located  between  the  P-  and 
S-fronts,  and  one  between  the  S-front  and  the  surface. 


If  the  surface  load  lies  below  the  limit  p  ,  but  above 


another  value  p£  <  p^  ,  so  that 


Pe  '  '  Pl 


(7) 


only  one  of  the  two  plastic  regions  occurs,  which  one  occurs 
depends  on  the  parameters  v  and  V/cp  .  If  these  parameters 
lie  in  Range  I  or  II  of  Fig.  A,  the  conf igurations of  Fig.  5 
or  6  apply,  respectively.  In  these  cases  two  of  the  neutral 
regions  shown  in  Fig.  3  coalesce  and  become  a  region  of 
uniform  stress  below  yield  as  noted  in  Figs.  5  and  6.  The 
demarkation  line  between  Ranges  I  and  II  in  Fig.  A  is  given 
by  the  relation 

2 


V/cn  = 


(1-2V) 


P  (1-V)(1-3V) 

If,  finally,  p  <  p  ,  the  solution  becomes  entirely  elastic, 

O  —  xi 

Fig.  7. 

One  detail  obtained  in  Ref.  [1]  for  the  step  pressure 
is  crucial  for  the  present  purpose.  If  p^  is  not  only  larger, 
but  very  much  larger  than  pT  ,  it  was  found  that  one  of  the 
two  plastic  regions  in  Fig.  3  is  of  much  greater  importance 
than  the  other.  This  means  that  the  energy  dissipated  and 
the  changes  in  stress  and  velocity  in  the  "important"  region 
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•  PARAMETERS  FOR  WHICH  TABLES  ARE  GIVEN 
IN  SECTION  IV. 
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range:  I  SOLUTIONS,  PE<P0< 
FIG.  5 


8 


CONFIGURATION  OF  ELASTIC  SOLUTIONS 
FIG.  7 
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are  very  much  larger  than  in  the  other  plastic  region.  In 
the  "important"  region  the  bulk  of  energy  dissipation  and  of 
the  changes  in  stress  and  velocity  occur  for  values  of  <p  very 
close  to  a  critical  angle  $  for  which  the  basic  differential 
equations  arc  capable  of  a  singular,  discontinuous  solution. 
The  critical  angle  is  connected  with  the  fact,  that  in  the 
von  Mises  material  considered,  plastic  shock  fronts  may  exist, 
having  a  velocity 


_2  2(l+v)  G 

C  =  3  (  1  —  2  V  )  p 


(9) 


As  discussed  in  Ref.  [1],  at  such  shock  fronts  the  state  of 
stress  ahead  and  behind  the  front  must  be  at  yield,  but  -  in 
addition  -  the  major  principal  stress  must  be  at  a  right  angle 
to  the  front  arid  the  other  two  principal  stresses  must  be 
equal.  It  is  shown  in  Ref.  [1]  that  at  such  a  plastic  shock 
front  all  three  principal  stresses  change  by  the  same  amount 
A a,  while  the  component  of  the  velocity  normal  to  the  front 
changes  as  required  by  conservation  of  momentum. 


In  the  present  steady-state  problem  the  angle  <J>  may  be 
obtained  from  the  value  for  c, 


r 

2 ( 1+V)  G 
3 (l-2v)  p 


(10) 


The  location  of  this  front  depends  on  the  value  of  Poisson's 
ratio  V.  If  v  >  1/8  the  critical  angle  $  lies  between  4>p 
and  4>  g  • 
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While  in  the  actual  nroblem  the  above  stated  conditions 
for  a  discontinuity  are  not  exactly  satisfied  and  no  shock 
front  occurs,  the  conditions  arc  nearly  satisfied  when 
>>  p,  As  a  result,  very  rapid  changes  were  found  in 
Ref.  [11  in  an  extremely  narrow  range  of  the  angle  <j>  near  $ 
when  pQ  exceeds  about  5.  (See  the  tables  in  Ref.  [2].) 

The  fact  that  for  large  values  of  p^  the  plastic  effects 
are  nearly  a  "shock  front"  at  <j>  =  <j>,  suggests  that  approximate 
solutions  could  be  obtained  by  compressing  (approximately) 
all  plastic  effects  into  a  discontinuity  in  normal  pressure 
and  velocity  at  $,  and  treating  the  material  elsewhere  (ap¬ 
proximately)  as  entirely  elastic,  even  if  the  yield  condition 
is  slightly  violated.  For  the  case  of  a  step  pressure  pQ  >  p^  , 
this  approximate  solution  for  V  >  1/8  would  have  the  con¬ 
figuration  shown  in  Fig.  8.  In  addition  to  the  elastic  dis¬ 
continuities  at  4> g  and  4>p  ,  there  is  an  assumed  discontinuity 
at  $ ,  consisting  of  a  change  of  all  three  principal  stresses 
by  an  amount  Aa  and  an  appropriate  change  of  normal  velocity. 

In  the  region  between  these  fronts  the  stresses  and  velocities 
are  uniform,  but  the  yield  condition  in  one  or  two  of  the 
regions  will  be  violated.  Such  violation  is  considered  ac¬ 
ceptable  provided  it  is  not  too  severe. 

In  view  of  the  fact  that  correct  solutions  for  the  step 
pressure  are  available,  there  is  no  need  to  consider  this 
case  further,  but  the  equivalent  approximations  ca.i  be  made, 
and  will  be  made  below  for  the  case  of  a  decaying  surface  load. 
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PLASTIC 

SHOCK 


CONFIGURATION 


FOR  APPROXIMATE  SOLUTION 

FIG.  8 
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Assuming  the  discontinuities  stated  above,  at  <j)  ,  <j>s  and  <j> , 

it  will  be  assumed  that  the  elastic  differential  equations 
apply  in  each  of  the  three  wedge  shaped  regions  shown  in 
Fig.  8.  Solving  the  boundary  value  problem  nosed  by  the 
prescribed  surface  load  and  by  momentum  or  continuity  con¬ 
siderations,  approximate  solutions  are  obtained.  After 
finding  the  solutions  the  degree  of  violation  of  the  yield 
condition  can  be  checked  by  computing  the  ratio, 


starting  near  the  leading  edge  of  the  load.  If  the  solution 
ir.  this  location  and  within  some  distance  r  from  the  wave 
front  is  acceptable,  the  hyperoolic  nature  of  the  problem 
permits  one  to  accept  the  solution  in  the  range  r  <  r,  even 
if  serious  violations  of  the  yield  condition  occur  elsewhere 
for  r  >  r.  A  further  restraint  on  the  range  of  validity,  r, 
is  the  fact  that  the  compressive  stress  at  the  shock  front 
<j>  =  $  must  inherently  increase  and  the  solution  can  only  be 
accepted  within  a  radius  r  where  this  requirement  holds. 

In  the  numerical  solution  obtained  later  the  ratios  n 
indicating  compliance  or  violation  of  the  yield  condition  are 
stated  for  each  point  where  stresses  are  given.  For  engineering 
purposes  values  up  to  n  =  1.2  were  considered,  but  a  reader  may 
elect  not  to  use  the  full  range  of  the  solution  given  if  smaller 
values  of  n  appear  appropriate.  It  is  found  that  the  range  of 
applicability  shrinks  as  V/Cp  becomes  smaller  and  approaches 
unity  . 
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The  formulation  of  the  problem  given  below  applies  only 
for  values  V  1/8.  Therefore  the  tables  given  contain  only 
values  V/Cp  j>  1.5  and  V  >  1/8. 
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SECTION  II 


(16) 


(17) 


3t  ,  y  3v 

37  +  IT  "  P  3t 


must  be  satisfied.  Equations  (12,  13,  15)  may  be  integrated 
with  respect  to  time,  requiring  the  addition  of  arbitrary 
functions  f^(x,y)  of  x  and  y, 


3u  1— V  * 


3x  E  x 


[0..  - 


fy  -  -lir  °x  +  V  + 


2  [7y  +  3x  2G  T  f3(x,y) 


while  Eq.  (14)  after  integration  may  be  written 


a  =  v(o  +  a  )  +  f  (x,y)  0 

z  x  y  o  w 

The  last  equation  is  the  only  relation  containing  0 z  ,  and 
simply  defines  this  quantity.  It  is  noted  that  the  arbitrary 
functions  ^  ,  f2  and  f3  are  not  entirely  "arbitrary",  but 
are  related  by  a  compatibility  relation,  as  may  be  seen  by 
eliminating  u  and  v  from  Eqs.  (18),  (19)  and  (20). 

The  solution  to  the  differential  equations  (16-20)  may 
be  expressed  in  terms  of  potentials  $  and  4*  which  satisfy  the 
wave  equations 


*  + 


—  $ 

?  tt 


4*  +4' 

xx  yy 


2  tt 


f 


Subscripts  on  potentials  indicate  derivatives.  In  terms  of 
these  potentials  and  of  two  arbitrary  functions  g  ^  and  g ^  of 
x  and  y,  the  displacements  become 


U  =  ^x  "  %  + 

v  =  *y  +  ^x  +  s2(x»y) 

while  the  velocities  u  and  v  are 

u  =  *  -  V 

x  y 

v  =  $  +  T 

y  x 


The  stresses  are 


o  =  2G  [7^-  <J>  +  ——  1-  -  T  ] 

x  1-2V  XX  1-2V  yy  xy 


O  -  2G  [-r~-  $  +  — 4>  +  Y  ] 

y  1-2V  xx  1-2V  yy  xy 


T  =  2G  [4  +  i  (Y  -  V  )] 

xy  2  xx  yy 


ith  a  given  by  Eq.  (21) 


( 2  A ) 

(25) 

(26) 

(27) 

(28) 

(29) 

(30) 


Substitution  of  Eqs.  (24—25)  and  (28-30)  into  Eqs. 
(18-20)  furnishes  relations  between  the  three  "arbitrary" 


functions  f^  ,  and  f^  ,  and  the  two  functions  g^  and  g ^  , 
i.e.  , 


fl  = 


3  x 


9g?  i  9g? 

- -  f  =  —  ( - i.  + - — 

3  y  ’  3  2  v  3y 


3  x 


) 


(31) 
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These  three  relations  and  the  Eqs.  (28-30)  satisfy  com¬ 
patibility  indentically,  leaving  the  two  functions  g  ^ (x , y) 
and  g2(x,y)  as  arbitrary  functions  in  lieu  of  the  related 


functions  £ ^  ,  £^  and  f^ 


It  is  now  noted  that  the  velocities  u  and  v,  and  the 
stresses  cry  ,  a  and  x  in  terms  of  the  potentials  4>  and  ¥ 
have  exactly  the  same  form  as  in  an  elastic  Diane  problem 
and  these  quantities  will  therefore  have  the  same  values  as 
in  an  entirely  elastic  problem  with  the  same  boundary  con¬ 
ditions.  The  plastic  deformations  at  the  shock  front  lead, 
however,  to  different  expressions  for  u,  v  and  0^  because  the 

functions  g,  ,  g„  and  f  now  occur.  As  will  be  seen  later, 
x  i.  o 

these  three  arbitrary  functions  follow  from  conditions  at  the 
P-front  and  at  the  plastic  front. 

Due  to  the  fact  that  the  steady-state  problem  is  con¬ 
sidered,  the  surface  pressure  is  solely  a  function  of 


£  -  x  -  Vt 

and  the  solutions  are  only  functions  of  £  and  y.  The  wave 
equations  for  the  potentials  then  become 


>yy  *  (M  '  X)  *es 


(32) 


(33) 


\y  -  <Ms  -  »  hi 


where 


M 


V_  .  V_ 

CP  5  s  "  c s 


18 


(34) 


I 


Noting  that  the  open  functions  f  ,  g  ^  and  g^  ,  which 


are  not  functions  of  t,  can  in  the  steady-state  not  depend 
on  x  either,  the  expressions  for  displacements  and  stresses 
become 


o  = 


a  = 


♦€  -  %  +  h(y) 

(35) 

*y  +  ^  +  82(y) 

(36) 

G[2(1  +  i?2„>  V  -  2V 

(37) 

G((M2  -  2)  *£5  +  2f£y] 

(38) 

G[2^y  -  <M2  -  2)  TE£] 

(39) 

v(a  +  a  )  +  f  (y) 
x  y  o 

(40) 

To  formulate  the  boundary  conditions,  expressions  will 
be  required  for  the  normal  and  tangential  velocities  u^  , 
and  the  normal,  tangential  and  shear  stresses  0^  ,  ,  T, 

for  planes  inclined  at  any  angle  $  (see  Fig.  1) 


N 

V 


_N 

G 


i  d>  <5> -  cos  4>  <J>  r  -  cos  <j>  T-p  -  sin  4>  (41) 

?y  C?  ?y 

-  cos  4>  -  sin  4>  4>_  -  sin  4>  H*,..  +  cos  4>  (42) 

55  5y  55  5y 


(Mg  (1  -  sin^)  -  2  cos  2<{>)  *  ^ 


-  2  sin  24.  4>?y  +  2  cos  24>  +  (M*  -  2)  sin  24»  ^  f  (43) 


5y 


55 
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2  1  _  o  v>  2 

=  [Mg  (1  -  y~-  cos  <J>)  +  2  cos  24>]  «J>^  + 

+  2  sin  2<J>  -  2  cos  2$  Y-  -  (M*  -  2) 

Cy  Cy  s 


sin  2<j>  V__ 
S 


(44) 


~  =  (M2  -  2)  sin  2 <f>  <3>r_  +  2  cos  24>  4>r  + 

r>  CC  Sy 


+  2  sin  2<j>  -  (Mg  -  2)  cos  2<fr 


(45) 


To  state  the  boundary  conditions  of  the  Droblem  properly, 
and  to  be  able  to  make  statements  on  the  character  of  the 
solutions,  the  differential  equations  for  <t>  and  ¥  must  be 
discussed.  While  both  differential  equations  are  hyperbolic, 
the  characteristics  for  $  are  steeper  than  the  boundary  ^ 
defining  the  plastic  front,  <[>p  <  $.  The  solution  for  $  in 
the  region  of  interest,  tji  <  <p  <  n,  Fig.  8,  is  therefore 

defined  by  statements  on  the  boundaries  <p  ~  <j>  and  <i>  =  ir . 

<I>  will  be  a  continuous  function  if  the  prescribed  conditions 
at  the  boundaries  are  continuous.  This  being  the  case,  it 
will  be  possible  to  use  for  $  an  expansion  in  a  Taylor  series 

in  £  and  y  for  a  region  near  the  origin  5  =  y  =  0. 

On  the  other  hand,  the  characteristic  directions  for  the 
differential  equation  on  ¥  are  flatter  than  the  boundary, 

<p  >  $.  The  differential  equation  does  therefore  not  apply 
at  <p  -  ,  and  the  solutions  V  for  4>  >  Vg  >  and  T  for 

$  <  <J>  <  4>„  are  not  fully  defined  by  prescriptions  on  <p  =  tt 

u 

and  <{)  =  <j> .  The  two  functions  T  and  ?  must,  however,  lead  to 
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values  of  stress  and  velocities  which  satisfy  momentum  con¬ 
siderations  at  4>  =  <f>  g  ,  which  furnish  additional  conditions 
for  their  determination.  While  at  least  some  of  the  deriva¬ 
tives  of  T  and  ?  at  <j>  =  <j>gWill  not  be  continuous,  each  of 
these  potentials  in  its  region  will  be  continuous,  and  each 
of  the  potentials  T  and  ?  can  again  be  expanded  in  a  Taylor 
series  . 

The  functions  $ ,  and  ?  are  therefore  subject  to  the 
conditions  listed  below.  On  the  surface  the  normal  and 
tangential  stresses  are  prescribed.  Thus,  for  an  exponentially 
decaying  surface  load,  the  stresses  are 

at  y  =  0 ,  £  <  0 


a  p_ 

a  _  1°  ^ 


(U  >  0) 


1,0 


(46) 

(47) 


At  the  plastic  front  <j>  =  $  conservation  of  momentum  requires 


at  <p  =  <t> 


AJ. 


4s 

k  N 


At 

k 


-  0 


A°N  =  AaT  =  A°z 


AO  =  j  AJ  x 


(48) 

(49) 

(50) 
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where  is  the  invariant 

J .  =  o  +  a  +  o  (51) 

1  x  y  2 

and  where  the  symbol  A  indicates  the  jump  in  stress,  or 
velocity  at  a  front  of  discontinuity  and  where  c  and  $  are 
given  by  Eqs.  (9)  and  (10),  respectively.  Moreover,  energy 
considerations  require  that  the  signs  of  Aa  and  of  the  total 
normal  stress  a N  at  $  are  equal. 

Due  to  the  fact  that  different  expressions  ?  and  ?  for 
the  solution  of  the  second  of  Eqs.  (33)  will  be  used,  it  is 
necessary  to  consider  conditions  at  the  location  <j>  =  <j><,  . 

Momentum  considerations  require  that  discontinuities  in  the 
shear  stress  and  the  tangential  component  of  the  velocity 
uT  are  proportional 

at  <j>  =  <t>s 

Atn  =  p  cg  Aut  (52) 

while  all  the  normal  and  tangential  stresses  and  the  normal 
velocity  at  (J>  *  <j)g  must  be  continuous 

at  <j>  =  <f>s 


Aon 

=  0 

(53) 

Aot 

=  0 

(54) 

Ac 

z 

=  0 

(55) 

aan 

=  0 

(56) 
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In  the  assumed  approximate  solution,  the  plnstic  shock 
propagates  into  a  region  which  has  been  previously  stressed 
by  the  passage  of  the  P-front.  In  this  region,  4>  p  <  <f>  <  <j> , 
the  stresses  are  uniform  and  just  at  the  yield  limit. 

9 

Their  values  are: 

for  <j>p  <  <j>  <  | 


ox 

-  /3  (1-v) 

1 1  —  2  sin^ijig  cot2^)p] 

(57) 

T 

l-2v 

a 

_X  _ 

-  /3  (1-v) 

cos  2({) s 

(38) 

k 

1-  2v 

T 

_2LZ 

k 

2/3  (1-v) 
1-  2v 

2 

sin  4> s  cot  <{>p 

(59) 

H*Q 

ii 

2  sin  2(j> 

P 

(60) 

The  above  equations  complete  the  available  conditions 
for  the  determination  of  t.ie  potentials  4> ,  *P  and  ? .  After 
the  potentials  are  obtained  the  arbitrary  function  f  (y ) 
appearing  in  Eq.  (40)  can  be  obtained  from  Eq.  (50).  The 
arbitrary  functions  g^(y)  and  §2^)  in  Eqs.  (35)  and  (36) 
are  defined  by  the  fact  that  displacements  at  the  plastic 
front  may  be  computed  from  the  velocities  in  the  region 
between  4> p  and  using  the  relations: 

for  <j>p  <  <S>  <  $ 

u  =  -  ^  (5  +  y  cot  <i>p)  (61) 

v  =  -  ^  (5  +  y  cot  <(>p)  (62) 
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SECTION  III 


SOLUTION  OF  THE  DIFFERENTIAL  EQUATIONS  BY  TAYLOR  SERIES 


The  differential  equations  for  the  potentials  to  be 
solved  are 


$ 

yy 

(M2  - 

1) 

-  0 

(63) 

? 

yy 

,  2 

1) 

\ 

(MS  " 

«  0 

'  (64) 

f  - 
yy 

(M s  - 

1) 

hi 

*  0 

The  boundary  conditions,  Eqs.  (46-47)  on  the  surface  give 


at  y  =  0,  (  <  0 

<Ms  -  2)  +  2hy  ‘  -  ^  <65> 

2f?y  -  <M2  -  2)  *K  -  0  (66) 

While  the  conditions,  Eqs.  (48-50)  at  the  plastic  front  give 
three  equations 


at  y  =  C  tan  <p 


|  )(H2  -  1)  +  cot  *  l*5y  +  »?el  +  - 


/l  .  1 

2^2 


1  /  1-V 
il-Tv 


CO  5 


/  4s 

vy 


(67) 
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I 


{  (M2  -  2)  sin  2$  +  2  cos  2$  +  2  sin  2$  9 


£y 


Cy 


-  (Mg  -  2)  cos  2|  ?  }  =  0 


(68) 


y  {  (M^  -  2)  cos  2$  -  2  sin  2$  <J>^  +  2  cos  2$  4* 

+  (Mg  -  2)  sin  2$  9^}  =  -  y  cos  2  ($  -  <j>p) 

Finally,  the  five  conditions,  Eqs.  (52-56)  at  the 
lead  to  only  one  condition  on  9  and  9 


S-f  ront 


(69) 


at  y  =  f  tan  (j>( 


<¥«  -  V  +  ta"  *s<?5y  -  V  -  0 


(70) 


The  surface  load  is  expanded  in  a  series  in  £ 


1*5  V  jl  .  m  .. 

P  e  =  p  )  —  y  £ 

r  n  1  n  ^  m  ^ 


1  .  m 


o  m! 

m=0 


(71) 


aoplicable  for  5  <  0,  while  the  potentials  $ ,  4*  and  9  are 
expanded  in  double  series 


G 

k 


51 


(72) 


The  differential  equations  and  the  conditions,  Eqs.  (46-56) 
are  homogeneous,  i.e.,  they  contain  consistently  only  second 
derivatives  of  <J>,  9  and  9  with  respect  to  £  and  y. 
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Substitution  of  Eqs.  (72)  into  Eqs.  (63),  (64)  and  (66-70) 
results  therefore  in  linear  equations  for  the  values  a,  , 
b.  ,  c,  which  couple  only  those  terms  with  indices 

*  J  1  >  Xs 

m  =  i+Jl-2.  The  potentials  can  therefore  be  conveniently 
written  as  a  series  of  polynomials.  Each  set  of  polynomials 
of  class  m  can  be  determined  independently  from  the  poly¬ 
nomials  for  a  different  value  of  m. 


G4> 

k 


g£ 

k 


G? 

k 


l  $(m) 

m  =  0 

m+2 

-  I 

i=0 

a  (m) 
i 

c1 

m+2  -  i 

y 

(73) 

00 

V  vp(m) 
m=0 

^  (m) 

m+2 

-  1 

1  =  0 

,  (m) 
bi 

c1 

m+2  -  i 

y 

(74) 

y  |(m) 

m  =  0 

|(.m) 

m+2 
-  1 

c  (m) 
i 

51 

m+2-i 

y 

(75) 

The  case  m  =  0  corresponds  to  the  step  load  and,  while 

quite  simple,  is  treated  differently,  see  Appendix  I.  For 

each  value  m  >.  1  there  are  (3m  +  9)  unknown  coefficients 

a^"1^,  c<*>.  Substituting  the  second  derivatives  of 

Eqs.  (73-75)  into  the  conditions  on  the  boundaries  of  the 

regions,  Eqs.  (65-70),  yields  six  equations  for  each  value 
* 

m  >_  1 ; 

m 

(Ms  '  2)(m  +  2)(m  +  U  am+2  +  2(ra  +  L)  Cm+1  =  “  ~TU  (?6) 

2 (m  +  i)  am+^  -  (M2  -  2)(m  +  2)(m  +  1)  cm+2  =  0  (77) 

m+2  , 

l  (b,  -  c.)  j  (tan  4>q)m+2_J  -  0  (78) 

4=1  -J  J  b 


The  superscripts  used  in  Eqs.  (73-75)  have  been  omitted 
for  simplicity. 
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(m  +  1 ) (m  +  2  )  [  ( M  2  -  1)  a  0  +  cot  $  b  _  ]  + 

m+2  m+2 

+  (m  +  l)(tan  4> ) m  [cot  <j>  a  ^  +  b ^ )  + 
m+1  . 

+  l  (m  +  3  -  j ) ( tan  $)J_i  { [ (m  +  2  -  j ) (H2  -  1)  + 
j=2 

+  (j  -  1)  cot2  <j>]  am+3_j  +  [ (m  +  2  -  j)  cot  $  + 

+  (j  -  1)  cot  $]  bm+3_j}  =  0  (79) 

(m  +  2)(m  +  l)[N5am+2  +  Ngbm+2^  +  +  l)(tan  <f>)m[Nga1  +  l^b.^]  + 

m+l  _  .  , 

+  l  (m  +  3  -  j  )  ( tan  <jj)J~  {  [  (m  +  2  -  j)  N  + 
j  =2  5 

+  (j  -  1)  cot  $  N6 ]  am+3_ j  +  [(m  +  2  -  j)  Ng  + 

+  (j  -  1)  cot  $  N?]  bm+3_j}  -  0  (80) 

(m  +  2 )  (m  +  l)[Ngam+2  +  NiQbm+2^  +  (m  +  D(tan  $)m[N7a^  -  N6b1]  + 
m+1 

+  l  (m  +  3  -  j )  (tan  $)',“i  {  [  (m  +  2  -  j  )  NQ  + 

J-2  9 

+  (j  -  1)  N 7  cot  $]  am+3_j  +  [ (m  +  2  -  j)  N1q  - 

-  (j  -  1)  N6  cot  $]  bm+3_1}  =  0  (81) 
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where 


N x  =  Mg  -  2  4-  sin2$  (4  -  2M2  -  M2) 

M2  -  4 

N2  =  -  (-^ - )  sin  2$ 

N3  -  2  +  (M2  -  4)  sin2f 
3M2  -  4 

N4  =  ( - ^ - )  sin  2$ 

N5  =  (M2  -  2)  sin  2$ 

N,.  -2  cos  2$ 

b 

N?  =  2  sin  2$ 

Ng  =  -  (Mg  -  2)  cos  2$ 

„r  .  ,..2  „2V  ,  1  +  2a/3  \  1  „2  ,  w2  ,  2r, 

Nn  =  2  [cos  2<J)  +  (M  -  M  )( - )  -  -r  M  +  M  sin  <{>] 

9  S  2  +  a/3  1 

N IQ  =  -  (M2  -  2)  sin  2$ 


(82) 


Substituting  Eqs.  (73-75)  into  the  wave  equations,  Eqs.  (63-64), 
2nd  s  a  u  3.  t  i  n  g  coefficients  of  like  powers  of  £  snd  y,  yields 
(3m  +  3)  recurrence  equations 
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i 


am-2q+2  =  <!<2  '  1  > "  am+2  ^  *„,l 


3m-  2q+l  '  <M  "  1>(I  am+l  62q-2,m  ^  f . ,  l 


bm-2q+2  =  (MS  '  l>’  V2  ^  P„,2 


b  =  f  m  ^  _  n'! 


m- 2q+l  S 


=  (m:  -  i)H  b  ,,  6 


q 

n  p 


m+1  2q-  2 ,  m  1  m ,  & 


o  q 

C  +  9  =  <Mc  “  Dq  C  ,0  n  P  # 

m-2q+2  S  m+2  m , 9, 


Cm-2q+l  (MS  “  1>q  Cm+1  62q-2,m  ^  Pm ,  9. 


\  (83) 


j 


where 


p  _  (m  +  3  -  21) (m  +  4  -  21) 
m ,  9,  21(29.  -  1) 


(m+2  -  29.)  (m  +  3  -  21) 
29.(29.  +  1) 


6 


2q-2 ,  m 


1 

0 


if  2q  -  2  i  m 
if  2q  -  2  =  m 


}  (84) 


i  <  q 


<  2m  +  3  +  (~l)m 
4 


j 


Equations  (76-81)  and  Eqs.  (83)  give  the  required  number  of 
(3m  +  9)  equations  for  the  (3m  +  9)  unknown  coefficients 
b^  and  c^  . 


a  . 

l 


Substitution  of  Eqs.  (83)  into  Eqs.  (76)  to  (81)  reduces 
the  system  to  a  set  of  six  linear,  nonhomogeneous ,  simultaneous 
equations,  involving  only  the  unknowns  am+2  >  am+l  ’  ^m+2  ’ 
bm+i  ,  cm+2  .  cm+i  •  After  solving  for  these  values,  the 
stresses,  velocities  and  displacements  for  each  of  the  m  com¬ 
ponents  of  the  surface  load  may  be  computed  at  any  location 
(£,y)  behind  the  front  at 
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SECTION  IV 


NUMERICAL  results  and  discussion 

Numerical  results  from  the  relations  derived  in  the 
previous  section  were  obtained  on  a  digital  computer, 
truncating  the  series  solution  by  retaining  terms  up  to 
and  including  m  =  7.  Retention  of  terms  m  >  7  would  have 
required  the  use  of  double  precision,  in  spite  of  the  fact 
that  the  computer  employed  carries  14  digits.  The  stresses 
found  are  presented  in  Tables  II  to  XXVIII  and  in  the  ac¬ 
companying  figures,  for  all  combinations  of  the  velocity 
ratios  V/cp  =  1.5,  2.0,  4.0  and  of  Poisson's  ratio  V  =  0.15, 
0.25,  0.35,  and  for  three  peak  values  of  surface  loading. 

The  values  of  pQ  /  k  are  selected  to  permit  interpolation  and 
extrapolation . 

The  approximation  employed  restricts  the  applicability 
of  the  approximate  solutions  to  distances  from  the  front  of 
the  applied  load  where  the  ratio  n,  Eq.  (11),  exceeds  unity 
only  moderately.  To  check  this  matter,  the  Tables  indicate 
the  value  n  for  all  points  where  stresses  are  listed. 
Truncation  of  the  series  further  limits  the  validity  of  the 
results  to  points  where  |u?|  <  1.7  and  only  such  points  are 
listed.  Finally,  the  range  of  applicability  of  the  present 
solution  is  limited  by  the  fact  that  the  approximate  results 
can  only  be  accepted  as  a  reasonable  approximation  in  the 
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range  where  the  shock  is  compressive.  The  tables  have  been 
restricted  in  range  such  that  this  requirement  is  never 
seriously  violated. 

For  each  set  of  values  V  and  V/Cp  considered,  results 

for  three  levels  of  the  load  pQ/k  are  given.  For  values  of 

p  /k  below  the  minimum  considered,  the  elastic  solution 
o 

presented  in  Appendix  II  may  be  used,  because  this  will 
generally  result  only  in  a  very  slight  violation  of  the  yield 
condition.  For  values  pQ/k  between  the  lowest  and  highest 
values  given  in  the  tables,  interpolation  between  the  three 
values  given  can  be  used  without  serious  error,  except  at 

•k 

point  (1  ).  For  values  pQ/k  which  exceed  the  highest  value, 

max(pQ/k)  considered,  the  principal  stresses  at  all  points, 

* 

except  at  point  (1  )  can  be  obtained  by  using  the  results 
for  the  respective  highest  value  of  pQ/k,  ana  adding  the 
increment 

0^in)  =  maxCp^/k)  -  PQ/k 

to  each  value  in  the  t-able  applicable  for  the  value 
max(pQ/k).  The  increment  is  inherently  negative, 

i.e.,  compressive.  Stresses  at  point  (1  )  are  always  those 
shown  in  the  tables  fqr  all  values  pQ/k  above  the  minimum 
considered . 

1 .  Discussion  of  typicaJ  results 

The  results  for  V/cp  =•  2.0,  v  =  0.25  for  the  three  loads 

p  /k  =  4.0,  6.0  and  10.0  are  discussed  as  typical, 
o 


(85) 
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The  principal  stresses  0j/k,  o^/k,  o  ^  /  k  ,  the  ans;le  0 
between  the  direction  of  and  the  C-axis,  and  the  values 
of  n  indicating  violation  of  the  yield  condition  are  calculated 
at  fourteen  points  on  a  rectangular  grid  in  the  £-y  plane. 

The  results  for  these  cases  are  given  in  Tables  XIV  -  XVI. 

In  every  table  point  1  has  three  sets  of  stresses  because 
the  stress  field  there  has  a  directional  singularity,  i.e., 
the  values  depend  on  the  direction  of  approach  <p ,  and  the 
results  differ  depending  on  whether  <j>p  <  4>  <  $ ,  $  <  <  <J>S  or 

<J>  <  <p  <  it.  The  three  sets  of  values  at  point  1  immediately 
behind  the  front  of  the  anplied  load  are  equal  to  those  due 
to  a  step  load  without  decay.  In  each  of  these  three  tables 
the  yield  violation  at  point  1  is  the  same,  n  =  1.056.  The 
yield  violations  for  a  step  wave,  for  a  range  of  the  para¬ 
meters  V/Cj.  ,  v  and  P0/k  are  shown  in  Fig.  9.  For  each  com¬ 
bination  of  V/Cp  and  V  there  is  a  limiting  load,  above  which 
the  degree  of  violation,  n,  remains  constant.  This  is  due  to 
the  fact  that  the  effect  of  any  higher  value  of  po/k  results 
only  in  a  jump  in  the  stress  invariant  at  <}>  =  <t> ,  a  matter 
which  does  not  affect  the  von  Misas  yield  condition. 

In  order  to  judge  the  seriousness  of  the  yield  violation 
at  point  1,  n  =  1.056,  Table  I  gives  the  exact  values  of  the 
stresses  O./k  and  of  the  angle  0,  for  a  load  of  p  /k  =  4.41, 

I  o 

as  obtained  in  Ref.  [1]. 
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TAELE  la  0  -  0.25,  V/c„-  2.0,  p  =  4. Aik 

P  o 

(Exact  Stresses) 


Locat ion 

-<Vk 

-  72/k 

-03/k 

8 

4>p  <  <P  <  <t>1 

2 .598 

0.866 

0.866 

60.00° 

CO 

•e- 

V 

-&• 

V 

CM 

-e- 

4.20 

2.47 

2.47 

67.46° 

4>s  <  <t>  <  $3 

4.20 

2.47 

2 .47 

78.98° 

<K  <  <l>  <  n 

4.41 

2.63 

2.73 

90.00° 

The  comparable  stresses  in  the  approximate  solution  for 

p  =  4.41k  are 
ro 


TABLE  lb 


V  =  0.25,  V / cp  =  2.0,  P( 
(Approximate  Stresses) 


=  4.41k 


Location 

-a^/k 

-°2/k 

-<r3/k 

0 

4>p  <  <l>  <  $ 

2.598 

0.866 

0.866 

Ov 

O 

o 

o 

o 

<t>  <  <j>  <  (j> 

4.345 

2.613 

2.613 

60.00° 

O 

4>s  <  <p  <  4> 

4.41 

2.549 

2.549 

90.00° 

_ _ _ 

The  differences  in  the  principal  stresses  in  the  two 
tables  are  not  large,  the  maximum  being  7%.  While  the 
values  8  listed  in  Table  la  agree  with  those  in  Table  lb 
it  is  noted  that  in  the  exact  case  the  value  8  varies 
between  (|>  <  <j>  <  <j> ^  from  60.00°  to  67.46°  and  between 

4> 0  <  <j>,  from  78.98°  to  90.00°.  In  the  approximate 


PRINCIPAL  STRESS  6  A  AT  4> 
FOR  V«  0.15  v/c„sl.50 

FIG.  10 


TABLE  II 


RESULTS  FOR  V  =  0.150  V/c„  =1.5  p  =  4.0k 

P  o 


^P 


Pt. 

y? 

yy 

O-j^/k 

Cfj/k 

°3/k 

o 

<D 

n 

1* 

0.000 

0.000 

-2.103 

-0.371 

-0.371 

48.189 

1  .uuu 

1** 

0.000 

G.UOO 

-3.72v, 

-1.987 

-1.987 

48.189 

1  •  JOO 

1*** 

0.000 

0.000 

-4.000 

-1.707 

-1.987 

90.000 

1.25u 

2 

-0.099 

0.000 

-3.619 

-1.613 

-1.916 

90.000 

1.081 

3 

-0.199 

0.000 

-5,27b 

-1.522 

-1.851 

90.000 

0.9^1 

4 

-0.299 

0.000 

-2.964 

-1 .434 

-1. /91 

90.000 

U.  800 

5 

-0.099 

0.050 

-3.578 

-1.846 

—  1 .  b46 

48.189 

1 .000 

6 

-0.199 

0.050 

-3.509 

-1.460 

-1.778 

89.645 

1.102 

7 

-0.299 

0.050 

-3.173 

-1.379 

-1.715 

89.748 

0 « 9  8^ 

8 

-0.199 

C.  100 

-3.443 

-1.710 

-1.710 

48.195 

1.000 

9 

-0.299 

0.100 

— 3 . 404 

-1.313 

-1.645 

8  9*311 

1.123 

1C 

-0.299 

0.150 

-3.313 

-1.580 

— 1 . 580 

48.195 

1 . 0  4  0 

1 l***# 

-0.199 

0.094 

-3.737 

-1.395 

-1.718 

89.^00 

1.269 

12**** 

-0.299 

0.  141 

-3.615 

-1.249 

-1.591 

88.846 

1.2  7s 

‘J  3 

-1.500 

0 » ouo 

-0.890 

-0.576 

-1.351 

90 , 0  0 \j 

u .  39o 

* 

For 

*P 

< 

<t>  < 

*  * 

For 

$ 

< 

V 

*#* 

For 

*S 

< 

4>  1 

TT 

**** 

At 

♦t 

At  *  1.357  k 
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TABLE  III 


RESULTS  FOR  v  a  0.150,  V/cp  *  1,5,  Pq  .  6.0k 


Pt. 

P? 

Py 

rP 

a1/k 

a2/k 

a^/k 

6° 

n 

1* 

0.000 

0.000 

-2.103 

-0.371 

-0.371 

48.189 

1.000 

1** 

0.000 

0.000 

-5.720 

-3.987 

-3.98? 

48.189 

1.000 

1#** 

0.000 

0.000 

-6.000 

-3.707 

-3.987 

90.000 

1.250 

2 

-0.059 

0.000 

-5.651 

-3.622 

-3.922 

90.000 

1.094 

3 

-0.1J9 

0.000 

-5.321 

-3.538 

-3.860 

90.000 

0.950 

4 

-0.179 

0.000 

-5.012 

-3.456 

-3.802 

90.000 

0.817 

5 

-r  059 

0.030 

-5.591 

-3.859 

-3.859 

48.189 

1.000 

6 

-0.119 

0.030 

-5.547 

-3.479 

-3.796 

89.630 

1.113 

7 

-0.179 

0.030 

-5.223 

-3.401 

-3.735 

89.668 

0.970 

8 

-0.119 

0.060 

-5.466 

-3.733 

-3.733 

48.199 

1.000 

9 

-0.179 

0.060 

-5.447 

-3.340 

-3.671 

89.278 

1.132 

10 

-0.179 

0.090 

-5.344 

-3.611 

-3.611 

48.198 

1.000 

11#**# 

-0.119 

0.056 

-5.759 

-3.421 

-3.741 

89.262 

1.267 

12**** 

-0.179 

0o085 

-5.644 

-3.284 

-3.622 

88.923 

1.276 

13 

-1.500 

0.000 

-2.010 

-1.335 

-3.033 

90.000 

0.855 

* 

For 

< 

<p 

< 

-©- 1 

M  M 

For 

£ 

T 

< 

< 

rf> 

'  s 

*** 

For 

*s 

< 

4> 

< 

IT 

#*** 

At 

*s 

Ax  -  1.357k 


TABLE  IV 


Pt . 

US 

uy 

0,  / k 

1* 

0.000 

0.000 

-2.103 

1** 

0.000 

0.000 

-9.720 

1*** 

0.000 

0.000 

-10.000 

2 

-0.039 

0.000 

-9.608 

3 

-0.080 

0.000 

-9.231 

4 

-0.119 

0.000 

-8.869 

5 

-0.039 

0.020 

-9.576 

6 

-0.080 

0.020 

-9.490 

7 

-0.119 

0.020 

-9.118 

8 

-0.080 

0.040 

-9.436 

9 

-0.119 

0.040 

-9.375 

10 

-0.119 

0.060 

-9.298 

11#*#* 

-0.080 

0.037 

-9.729 

12**#* 

-C.119 

0.056 

-9.599 

o2/k  a3/k  0° 


-0.371 

-0.371 

48.189 

1.000 

-7.987 

-7.987 

48.189 

1.000 

-7.707 

-7.987 

90.000 

1.250 

-7.612 

-7.914 

90.000 

1.075 

-7.518 

-7.844 

90.000 

0.909 

-7.425 

-7.776 

90.000 

0.753 

-7 . 844 

-7.844 

48.189 

1.000 

-7.451 

-7.772 

89.552 

1.096 

-7.361 

-7.703 

89.550 

0.931 

-7.702 

-7.702 

48.205 

1.001 

-7.292 

-7.632 

89.126 

1.117 

-7.564 

-7.564 

48.205 

1.000 

-7.387 

-7.711 

89.173 

1.269 

-7*230 

-7.576 

88.787 

1.279 

*  For<j>p<4><4> 

##  For  $  <  <j>  <  <PQ 

J 

*##  For  4>  <  <_  ir 

»  ■  + 

*#**  At  4> 

O 

At  =  1.357k 
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PRINCIPAL  STRESS  6%  AT  ? 
FOR  V  *  0.15  V/Cp*2.0 


h-  y 

tan^» 


TABLE  V 

RESULTS  FOR  v  =  0.150,  V/cp  =  2.0,  p 


3.0k 


4  3  2  1 


rP 


Pt.  p? 


1* 

0.000 

1** 

0.000 

1*** 

o.oou 

2 

“0.281 

3 

-0.561 

4 

-0.842 

c 

“0.281 

6 

“0.561 

7 

-0.842 

8 

-0.561 

V 

-0.842 

10 

-0.842 

11**** 

-0.561 

12**** 

“0.842 

* 

For 

*P 

< 

** 

For 

$ 

< 

<t> 

*** 

For 

*S 

< 

4> 

**** 

At 

♦t 

uy 

Oj/k 

0.000 

“2.103 

0.000 

-2.859 

0,000 

-3.000 

0.000 

— 2 . 266 

0.000 

“1.712 

0.000 

-1.293 

0.100 

-2.600 

0.100 

-2.048 

0.100 

-1.540 

0.200 

-2.367 

0.200 

“1.855 

0.300 

“2.157 

0.190 

“2.430 

0.285 

“2«  190 

< 

-e-i 

< 

<i  _ 

b 

< 

n 

o2/k 

o3/k 

-0.371 

“0.371 

“1.127 

“1.127 

-0,986 

-1.127 

-0.783 

“0.986 

-0.605 

-0.877 

-0.448 

-0.790 

-0.868 

-0.868 

-0.576 

-0.742 

-0.425 

-0*643 

-0.635 

“0,635 

-0.387 

-0.521 

-0.425 

-0.425 

-0.531 

—0 . 644 

-0.337 

-0.438 

6° 

n 

60.000 

1.000 

60.000 

1.000 

60.000 

1.124 

90.000 

0,804 

90.000 

0.57  / 

90.000 

0.425 

60.000 

1.000 

91.690 

0.806 

92.970 

0.591 

60.000 

1  .OoO 

93.180 

0.812 

60.000 

1.000 

91.330 

1.065 

92*000 

1.042 

At 


0.944k 


TABLE  VI 


RESULTS  FOR  V  «  0.160,  V/c. 


2.0  , 


6.0k 


NOT  TO  SC 


Pt. 

y? 

yy 

c^/k 

a2/k 

1* 

0.000 

0.000 

-2.103 

-0.371 

1** 

0.000 

0.000 

-5.859 

-4.127 

1**# 

0.000 

0.000 

-6.000 

-3.986 

2 

-1.122 

0.000 

-1.953 

-2.630 

3 

-1.403 

0.000 

-1.474 

-2.377 

4 

-1.683 

0.000 

-1.107 

-2.148 

5 

-1.122 

0.400 

-4.078 

-2.345 

6 

-1.403 

0.400 

-3.118 

-2.071 

7 

-1.683 

0.400 

-2.516 

-1.707 

8 

-1.403 

0.500 

-3.736 

-2.005 

9 

-1.683 

0.500 

-2.865 

-1.766 

10 

—1.683 

0.600 

-3.430 

-1.698 

n**** 

-1.403 

0.475 

-3.584 

-2.010 

12**** 

-1.683 

0.570 

-3.249 

-1.720 

03/k 


-0.371 

60.000 

1.000 

-4.127 

60.000 

1.000 

-4.127 

90.000 

1.124 

-3.317 

90.000 

0.682 

~3 • 20  » 

90.000 

0.867 

-3.117 

90.000 

1  006 

-2.345 

60.000 

1.000 

-2.160 

1U6.660 

0.580 

—2 . 016 

127. 04L 

U.4U8 

-2.005 

60.000 

1.000 

-1.838 

108.970 

0.615 

-1.698 

60.000 

1.000 

-2*040 

97.520 

0.900 

-1.736 

99.080 

0.879 

<({:<$ 

<  <j>  <  <j), 

v 

<  4>  <  tr 


TABLE  VII 


RESULTS  FOR  v  =  0.150,  V/cp  a  2.0,  pq  =  10.0k 


rP 


Pc  . 

yy 

Ox/k 

c?  2  /  k 

03/k 

0° 

n 

1* 

0.000 

0.000 

-2.103 

-0.371 

-0.371 

60.000 

1.000 

1** 

0.000 

0.000 

-9.859 

-8.127 

-8.127 

60.000 

1.000 

1#** 

0.000 

0.000 

-10.000 

-7.986 

-8.127 

90.000 

1.124 

2 

-0.280 

0.000 

-7,553 

-7,308 

-7.658 

90.000 

0.179 

3 

-0.561 

0.000 

-6,716 

-5.705 

-7.292 

90.000 

0.803 

4 

-0.701 

0.000 

-6.447 

-4.958 

-7.140 

90.000 

1.114 

5 

-0.280 

0.100 

-8.996 

-7.264 

-7.264 

60.000 

1,000 

6 

-0.561 

0.100 

-6.899 

-6.548 

-6.842 

117.688 

0.188 

7 

-0.701 

0.100 

-6.425 

-5.856 

-6.667 

71.592 

0.416 

8 

-0.561 

0.200 

“8.218 

-6.487 

-6.487 

59.989 

0.999 

9 

-0.701 

0.200 

-7.183 

-6.177 

-6.285 

102.388 

0.552 

10 

-On  701 

0.250 

-7.859 

-6.128 

-6.128 

59.989 

0.999 

1 1**** 

-0.561 

0.190 

-8.108 

—6 .460 

-6.518 

95.150 

0.935 

12**** 

-0.701 

0.237 

-7.703 

-6.122 

-6  *165 

96.651 

0.900 

13**** 

-1.500 

0.508 

-5.824 

-4 .466 

-4.466 

106.026 

0.784 

14 

-1.500 

0.000 

-5.168 

-2.225 

-6.538 

90.000 

2.203 

* 

For 

< 

<p 

< 

** 

For 

i-e- 

< 

<p 

< 

CO 

•e- 

*** 

For 

*s 

< 

<t> 

< 

IT 

**** 

At 

+s 

At  «  0.944k 
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TAb_E  VIII 


RESULTS  FOR  v  =  0.150,  V/cp  >  4.0.  PQ  •  3.0^ 


Ft  . 

V5 

py 

a  1  /  k 

02/k 

a3/k 

0° 

1* 

0.000 

0.000 

-2.103 

-0.371 

-0.371 

75.522 

0.000 

0.000 

-2.968 

-1.236 

-1.236 

75.522 

1*** 

0.000 

0.000 

-3.000 

-1.204 

-1.236 

90.000 

2 

-0.587 

0.000 

- 1,667 

-0.783 

-0.973 

90.000 

3 

-0.798 

0.000 

-1.349 

-0.661 

-0.907 

90.000 

4 

-1.174 

0.000 

-0.926 

-0.475 

-0.815 

90.000 

5 

-0.587 

0.100 

-2.493 

-0.761 

-0.761 

75.522 

6 

-0.798 

0.100 

-2.046 

-0.600 

—0 .670 

107.294 

7 

-1.174 

0.100 

-1.471 

-0.343 

-0.545 

113.343 

8 

-0.798 

0.136 

-2.342 

-0.612 

-0.612 

75.549 

9 

-1.174 

0.136 

-1.724 

-0.294 

-0.472 

114.592 

10 

-1.1/4 

0.200 

-2.102 

-0.372 

-0.372 

75.546 

11###* 

-0.798 

0.129 

-2.314 

-0.584 

-0.620 

108.183 

12**** 

-1.174 

0.190 

-2.146 

-0.255 

-0.383 

114.461 

13 

-1,500 

0.000 

—0 . 667 

-0.336 

-0.756 

90.000 

# 

For 

< 

<P 

< 

<j> 

** 

For 

-O-l 

< 

<P 

< 

*s 

*## 

For 

*S 

< 

<J> 

< 

n 

**** 

At 

♦t 

45 


n 


1.000 

1.000 

1.027 

0.465 

0.348 

0.235 

1.000 

0.815 

0.601 

0.998 

0.779 

0.999 

0.988 

1,056 

0,221 


At  =  0.442k 


NOT  TO  SCALE 


<j>p  *  165.52 
d>„  ■  170.77 


170.33 


Pt. 

u5 

vy 

c^/k 

a2/k 

a3/k 

0° 

n 

1* 

0.000 

0.000 

-2.103 

-0.371 

-0.371 

75.520 

1.000 

1*# 

0.000 

0.000 

-5.968 

-4.236 

-4.236 

75.520 

1.000 

1*** 

0.000 

o.ooo 

-6.000 

-4.204 

-4.236 

5'  .000 

1.028 

2 

-0.294 

0  e  000 

-4.4?3 

-3.747 

-3.939 

90.000 

0.377 

3 

-0.587 

0.000 

-3.361 

-3.335 

-3.710 

90.000 

0.209 

4 

-0.881 

0.000 

-3.302 

-2.486 

-3.533 

90.000 

U.524 

5 

-0.294 

0.050 

-5.470 

-3.738 

-3.738 

75.520 

1.000 

6 

-0.587 

0.050 

-4.214 

-3.196 

•'3.468 

119.020 

0.527 

7 

-0.881 

0.050 

-3.451 

-2.578 

-3.261 

49.700 

0.459 

8 

-0.587 

0.100 

-5.016 

-3.288 

-3.288 

75.520 

1.000 

9 

-0.881 

0.100 

-4.181 

-2.495 

-3.044 

128.680 

0.860 

10 

-0.881 

0.150 

^4.607 

-2.879 

-2.879 

75.520 

1.000 

11#*#* 

-0.587 

0.095 

-5.035 

-3.129 

-3.300 

116.430 

1.083 

12**** 

-0.881 

0.143 

-4.891 

-2.473 

-2.896 

123.800 

1.291 

13 

-1.500 

0.000 

-1.335 

-2.469 

-3.276 

90.000 

0.975 

For  <j>p  <  4>  <  <P 
For  |  <  (p  <  <f>, 


For  4>s  <  4>  £  tr 


###* 


At  <{), 


At  c  0.442k 


‘j  / 


TABLE  X 

RESULTS  FOR  v  =  U.150  V/c 


NOT  TO  SCALE 


<J>p  a  165.52 
4>s  =  170.77 
$  =  170. S3 


Pt . 

uK 

py 

0]L/k 

02/k 

a3/k 

0° 

n 

1* 

O.OOo 

0.  ouo 

-2.103 

-0.371 

-0.371 

75.520 

1 .  OOO 

1** 

O.OGo 

C  »  OoO 

-9.968 

-8.236 

-8.236 

90.000 

1  .  JuU 

1*** 

0.000 

0.OU0 

-lU.OOo 

-8.204 

-8.236 

90.000 

1  .U2B 

2 

-0.117 

0.000 

-0 . 892 

-7.883 

-8.022 

90.000 

u .  547 

3 

-0.235 

0.000 

-  / . 906 

-7.584 

-7.829 

9  0  s  0  0  0 

c  .  16^ 

4 

-0.352 

O.ooo 

-7 .305 

—7 .030 

-7.656 

9o  .woo 

w.3l‘- 

5 

-0.117 

0.020 

-9.626 

-7.894 

-7.894 

75.520 

1 . 0  0  o 

6 

-0.235 

0.020 

-8.586 

-7.558 

-7.688 

107.500 

u.560 

7 

-0.352 

C.020 

-7.765 

-7.145 

-7.503 

126. 46u 

0.311 

3 

-0.233 

0.040 

“9.292 

-7.568 

—  7 .568 

75.520 

1.000 

9 

-0.352 

0.040 

-8.441 

-7.099 

—  7.370 

119.720 

o.7lo 

10 

-0.352 

0.060 

-8.976 

-7.252 

-7.252 

75.520 

1  .OuO 

-0.235 

0.038 

-9.269 

-7.523 

-7.575 

109.480 

0.99  3 

12***# 

-0.352 

0.057 

-9.058 

-7.072 

-7.264 

117.250 

1.095 

13 

-0.650 

0.000 

-6.675 

-5.220 

-7.290 

90.000 

1.063 

* 

For 

< 

<P 

< 

5 

** 

For 

$ 

< 

< 

* 

*** 

For 

*S 

< 

0 

< 

TT 

**** 

At 

<j,+ 

At*  0.442k 


tan$> 


6j,  - 

PRiNCiPAL  SiRfcSS  'K  AT  <P 
FOR  "7=0.25  %p  =  1.5 

FIG.  13 
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TABLE  XI 


RESULTS  FOR  v  =  0.250,  V/cp  3  1»5,  pQ  »  ^.Ok 


1* 

1** 

1*** 

2 

3 

4 

5 

6 

7 

8 
o 

10 

1 1*### 
12**** 
13**** 

14 

15 


0.000 

0.000 

0.000 

-0.174 

-0.349 

-0.524 

-0.174 

-0.349 

-0.524 

-0.349 

-0.524 

-0.524 

-0.349 

-0.524 

-1.500 

-1.500 

-1.50C 


0-000 
0.000 
0.000 
0.000 
0.000 
0.000 
0.100 
0.100 
0.100 
0.200 
0.200 
0.300 
Oo  145 
0.218 
0.625 
Ok  000 
0.858 


-2.598 

-3.900 

-4.000 

-3.359 

-2.820 

-2.368 

-3.677 

-3.243 

-2.717 

-3.467 

-3.149 

-3.270 

-3.468 

-3.239 

-2.281 

-0.890 

-2.369 


-0.866 

-2.168 

-2.069 

-1.802 

-1.667 

-1.360 

-1.945 

-1.526 

-1.326 

-1.735 

-1.252 

-1.538 

-1.491 

-1.233 

-0.107 

-0.543 

-0.637 


— 0 . 866 
-2.168 
-2*168 
-1.941 
-1.748 
-1.583 
-1.945 
-1.732 
-1.550 
-1.735 
-1.535 
-1.538 
-1.731 
-1.534 
-0.652 
-1.009 
-0.637 


48.189 

48.189 

90.000 

90.000 

90.000 

90.000 

48.189 

87.788 

88.079 

48.189 

85.860 

48.189 

86.718 

85.430 

81.240 

90.000 

48.189 


* 

For 

< 

<t> 

< 

<J> 

** 

For 

$ 

< 

4* 

< 

♦s 

*** 

For 

*8 

< 

< 

TT 

**** 

At 

n 


1.000 

1.000 

1.087 

0*861 

0.677 

0.529 

1.000 

0.937 

0.746 

1.000 

1.023 

1.000 

1.079 

1.082 

1.131 

0.242 

1.000 


At  =  1.223k 


TABLE  XII 


RESULTS  FOR  V  *  0.250,  V/cp  =  1.5,  pQ  =  6.0k 


Pt. 

U5 

uy 

Oj  k 

C2/k 

CJ3/k 

0° 

n 

1* 

0.000 

0.000 

“2.598 

-0.866 

-0.866 

48.189 

1.000 

1** 

0.000 

0.000 

“5.900 

-4.168 

-4.168 

48.189 

1.000 

1#** 

0.000 

0.000 

-6.000 

-4.069 

-4. 168 

90.000 

1.087 

2 

-0.174 

0.000 

-5.038 

-3.669 

-3.828 

90.000 

0.748 

3 

“0.349 

0.000 

-4.230 

-3.316 

-3.538 

90.000 

0.476 

4 

-0.524 

0.000 

-3.552 

-3.005 

-3.290 

90.000 

0.273 

5 

-0.174 

0.100 

-5.565 

-3.833 

-3.833 

48.189 

1.000 

6 

-0.349 

0.100 

—4  it  867 

-3.251 

-3.513 

86.467 

0.867 

7 

“0.524 

0.100 

-4.077 

-2.952 

-3.241 

86.432 

0.584 

R 

“0. 349 

0.200 

-5.251 

-3.518 

-3.518 

48.189 

1.000 

9 

-0.524 

0.200 

-4.731 

-2.836 

-3.218 

83.742 

1.002 

10 

-0.524 

0.300 

-4.955 

-3.223 

-3.223 

48.189 

1.000 

11*#** 

-0.349 

0.145 

-5.208 

-3.198 

-3.512 

85.145 

1.U81 

12 **** 

-0.524 

0.218 

-4.868 

-2.806 

-3.217 

83.299 

1.091 

13#**# 

-1.200 

0.500 

-3.817 

-1.542 

-2.248 

78.788 

1.164 

14 

-1.200 

0.000 

-2.085 

-1.806 

-2.624 

90.000 

0.415 

* 

For 

-e- 

< 

<t> 

< 

-6-1 

*# 

For 

■©■I 

< 

<t> 

< 

##* 

For 

< 

<t> 

< 

U 

****  At 

O 


At  =  1.223k 
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NOT  TO  SCALE 


<{>p  =  138.18 
«j>s  =  137.36 
$  =  150.20 


Pt. 

y? 

yy 

a1/k 

02/k 

a3/k 

0° 

n 

1* 

0.000 

0.000 

-2.598 

—0 . 866 

-0.866 

48.189 

1.000 

1** 

0.000 

0.000 

-9.900 

-8.168 

-8.168 

48.189 

1.000 

1*#* 

0.000 

0.000 

-10.000 

-8.069 

-8.168 

90.000 

1.08  7 

2 

-0.174 

0.000 

-8.397 

-7.402 

-7.601 

90.000 

0.52b 

3 

-0.349 

0.000 

-7.051 

-6.814 

-7.117 

90.000 

0.159 

4 

-0.524 

0.000 

-6.295 

-5.921 

-6.705 

90.000 

0.392 

5 

-0.174 

0.100 

-9.342 

-7.610 

-7.610 

48.189 

1.000 

6 

-0.349 

0.100 

-8.122 

-6.697 

-7.077 

83.278 

0.737 

7 

-0.524 

0.100 

-6.811 

-6.192 

-6.623 

78.957 

0.317 

8 

-0.349 

0.200 

-8.817 

-7.085 

-7.085 

48.189 

l.OOu 

9 

-0.524 

0.200 

-7.910 

-5.989 

-6.584 

79.569 

0.9e3 

10 

-0.524 

0.300 

-8.325 

-6.592 

-6.592 

48.189 

1.000 

11***# 

-0.349 

0.145 

-8.694 

-6.602 

-7.074 

82.168 

1.097 

12*#** 

-0.524 

0.218 

-8.142 

-5.936 

-6.583 

79.413 

1.134 

13 

-1.060 

0.000 

-5.029 

-3.463 

-5.774 

90.000 

1*179 

* 

For 

< 

<t> 

< 

* 

## 

For 

i-e- 

< 

♦ 

< 

* 

##* 

For 

*S 

< 

* 

< 

IT 

****  At  <{>, 


At  «  1.223k 


SURFACE 


V 


i  ABLE  XIV 

RESULTS  FOR  V  =  0.250,  V/cp  «  2.0,  PQ  *  4.0k 


Pt . 

M? 

My 

1* 

0.000 

0.000 

1#* 

0.000 

0.000 

1#*# 

0.000 

0.000 

2 

-0.996 

0.000 

3 

-1.245 

0.000 

4 

-1.494 

0.000 

5 

-0.996 

0.400 

6 

-1.245 

0.400 

7 

-1.494 

0.400 

8 

-1.245 

0,500 

9 

-1.494 

0.500 

10 

-1.494 

0.600 

1  i###'* 

-1.245 

0.375 

12**** 

-1.494 

0.450 

13##*# 

-i.700 

0.513 

14 

-1.700 

0.000 

* 

For 

< 

< 

** 

For 

$ 

< 

$ 

< 

*s 

*#* 

For 

*s 

< 

4» 

< 

7T 

##*# 

At 

4>+ 
f  s 

At 

n  0 

• 

»o 

o 

o 

C^/k 

02/k 

-2.598 

-0.866 

-3.936 

2.203 

-4.000 

-2.139 

-1.477 

-1.026 

-1.151 

-0.847 

-0.896 

-0.693 

-2.960 

-1.228 

-2.502 

-0.769 

-1.871 

-0.674 

-2.758 

-1.026 

-2 . 33^ 

-0.603 

-2.572 

-0.840 

-2.329 

-0.831 

-2.093 

-0.640 

-1.916 

-0.496 

-0.725 

-0.579 

a3/k 

9° 

-0.866 

60.000 

-2.203 

60.000 

-2.203 

90.000 

-1.294 

90.000 

-1.168 

90.000 

-1.066 

90  o000 

-1.228 

60.000 

-0.999 

54.610 

-0.817 

89.640 

-1.026 

60.000 

-0.814 

55.000 

-0.840 

60.000 

-0.997 

87.040 

-0.812 

86.950 

-0o671 

86.990 

-0.995 

90.000 

3.000 

1.000 

1.056 

0.227 

0.181 

0.187 

1.000 

0.941 

0.654 

1.000 

0.944 

1.000 

0.821 

0.794 

0.774 

0.211 
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f 


I  RESULTS  FOR 


TABLE  XV 
v  =  0.230.  V/cp 


2*0,  po 


6.0k 


PC. 

P? 

Py 

<Vk 

a2/k 

a3/k 

6° 

n 

1* 

0,000 

0.000 

-2.598 

-0.866 

-0,866 

60,000 

1.000 

1*# 

0.000 

0.000 

-5.935 

-4*203 

-4.203 

60.000 

1.000 

1*** 

0.000 

0.000 

-6.000 

-4.139 

-4.203 

90.000 

1.056 

2 

-0.996 

0.000 

-2.469 

-2.216 

-2.840 

90.000 

0.313 

3 

-1.245 

0.000 

-2.201 

-1.726 

-2.650 

90.000 

0.462 

4 

-1.494 

0.000 

-1.969 

-1.343 

-2.497 

90.000 

0.577 

5 

-0.996 

0.400 

-4.471 

-2.739 

-2.739 

60.000 

1.000 

6 

-1,245 

0.400 

-3.796 

-2.040 

-2.396 

51.966 

0.928 

7 

-1.494 

0.400 

-2.806 

-1.941 

-2.123 

89.261 

0.456 

8 

-1.245 

0.500 

-4.169 

-2.437 

-2.437 

59.996 

1.000 

9 

-1.494 

0.500 

-3.542 

"1.792 

-2.119 

52.540 

0.930 

10 

-1.494 

0.600 

-3.890 

-2.158 

-2.158 

59.996 

1.000 

11*##* 

-1.245 

0.375 

-3.497 

-2.172 

-2.393 

84.958 

0./09 

12**** 

-1.494 

0.450 

-3.143 

-1.886 

-2.116 

84.693 

0 . 6  o  9 

13**** 

-1.700 

0.512 

-2.878 

-1.670 

-1.905 

84.666 

0 . 640 

14 

-1.700 

0.000 

-1.799 

-1.087 

-2.390 

90.000 

0.652 

* 

F-or 

< 

< 

$ 

** 

For 

i-e- 

< 

<P 

< 

*s 

###  For 

< 

* 

< 

TT 

**** At 

<j,+ 

At  *  0.900k 
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TABLE  XVI 


RESULTS  FOR  V  =  0.250,  V/cp  =  2.0,  pQ  =  10.0k 


Pt . 

0? 

py 

°l^k 

°2/k 

03/k 

0° 

n 

1* 

0.000 

0.000 

-2.598 

-0 .866 

-0.866 

60.000 

l.oou 

1** 

0.000 

0.0O0 

-9.935 

-8.203 

-8.203 

60 .000 

l.ouu 

1*** 

0.000 

0.000 

-10.UUU 

-8.139 

-8.203 

90.000 

1.056 

2 

-0.249 

0.000 

-7.795 

-7.241 

-7.428 

yu.uo-j 

0.281 

3 

-0.498 

0.000 

-6.499 

-6.076 

-6.812 

90.000 

0.369 

4 

-0.747 

0.000 

-5.879 

-4.737 

-6 .322 

9U.U0U 

0.817 

5 

-0.249 

0.100 

-y.254 

-7.522 

-7.322 

60.000 

l.ouu 

6 

-0.498 

0.100 

-7.301 

-6.599 

-6.803 

86.970 

0.361 

7 

-0.747 

0.  100 

-3.959 

—5.661 

—6 .233 

84.499 

0.266 

8 

-0.498 

0.200 

—8.622 

-6.890 

-6.890 

59.992 

1.000 

9 

-0.747 

0.200 

-6.869 

-3.979 

-6.224 

83.636 

G.40y 

10 

-0.747 

0.300 

-8.037 

-6.305 

-6.303 

59.993 

1.000 

1 1**** 

-0.498 

0.150 

-8.048 

-6.622 

-6.834 

84.925 

0.769 

12**** 

-0.747 

0.225 

-7.231 

-3 .967 

-6.236 

82.589 

0.666 

13**** 

-1.000 

0.3U1 

-6.491 

“3*365 

-5.680 

80.473 

0.380 

14 

-1.000 

0.000 

-3.347 

-3.678 

-5.925 

90.000 

1.166 

*  rOV<j)p<41''}' 

**  For  $  <  <J>  <  <}>s 

***  For  <  4>  <_  it 

b 

****  At  (jig 

At  *  0.900k 
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TABLE  XVII 


Pt . 

y5 

yy 

1* 

0.000 

0.000 

1** 

0.000 

0.000 

1*** 

0.000 

0.000 

2 

-0.527 

0.000 

3 

-1.054 

0.000 

4 

-1.581 

0.000 

5 

-0.527 

0.100 

6 

-1.054 

0.100 

7 

-1.581 

0.100 

8 

-1.054 

0.2C0 

9 

-1.581 

0.200 

10 

-1.581 

0.300 

11*#** 

-1.054 

0.153 

12**** 

-1.581 

0.230 

13#*#* 

-1.750 

0.255 

14 

-1.750 

0.000 

15 

-1.750 

0.331 

o1/k  C2/k  0; 

-2.598  -0.866  -0 
-3.983  -2.251  “2 

-4.000  -2.234  -2 

-2*360  -1.546  -1 

-1.393  -1.085  -1 

-0.819  -0.761  -1 

-3.485  -1.753  -1 

-2.112  -1.105  -1 

“1.379  -0.599  -0 

-3.053  “1.322  -1 

-2.027  -0.594  -0 

-2.678  -0.946  -0 

-2.604  -1.174  -I 

“2.244  -0.640  -0 

-2.159  -0.476  -0 

-0.687  -0.676  -1 

-2.569  -0.838  -0 


/  k  6° 


.866  75.522  1.000 
.251  75.522  1.000 
.251  90.000  1.014 
.669  90.000  0.439 
.312  90.000  0.159 
,087  90.000  0.174 
.753  75.522  1.000 
.248  111.589  0.545 
.938  124. 33y  0.390 
.322  75.524  0.999 
.883  121.034  0.  / 57 
.946  75.524  0.999 
.268  110.672  0.799 
.888  119.195  0.863 
.780  121.233  0.897 
.033  90.002  0.203 
.838  75.524  0.999 


* 

For 

Tl 

< 

$ 

< 

<P 

*# 

For 

I 

< 

<P 

< 

*s 

### 

For 

*S 

< 

<P 

< 

7T 

**## 

At 

♦t 

At  =  0.043k 
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TABLE  XVIII 


Pt. 

P? 

py 

Vk 

1* 

0.000 

0.000 

-2.598 

1** 

0.000 

0.000 

-5.983 

1*** 

0.000 

0.000 

-6.000 

2 

-0.527 

0.000 

-3.541 

3 

-1.054 

0.000 

-2.510 

4 

-1.581 

0.000 

-2.024 

5 

-0.527 

0.100 

-5.237 

6 

-1.054 

0.100 

-3.383 

7 

-1.581 

0.100 

-2,516 

8 

-1.054 

0.200 

-4,588 

9 

-1.581 

0.200 

-3.358 

10 

-1.581 

0.300 

-4.026 

11#*** 

-1.054 

0.153 

-4.071 

12**** 

-1.581 

0.230 

-3.646 

13**** 

-1.700 

0.248 

-3.575 

14 

-1.700 

0.000 

-1.934 

* 

For 

< 

<P 

< 

$ 

** 

For 

< 

<P 

< 

*#* 

For 

*S 

< 

<*> 

< 

IT 

**** 

At 

4>+ 

At  *  0.043k 


a2/k 

a3/k 

0° 

D 

'0.866 

-0.866 

75.522 

1.000 

■4.251 

-4.251 

75.522 

1.000 

'4.234 

-4.251 

90.000 

1.014 

■3.201 

-3.378 

90.000 

0,169 

■2.089 

-2.842 

89.999 

0.377 

1.225 

-2.506 

90.000 

0.645 

3.505 

-3.505 

75.522 

1.000 

2.326 

“2.747 

129.020 

0.532 

1.334 

-2.282 

56.390 

0.625 

2.857 

-2.857 

75.526 

0.999 

1.455 

“2.199 

133.128 

0.959 

2.294 

-2.294 

75.525 

0.999 

2.4?9 

-2.777 

121.415 

0.846 

1.563 

-2.207 

129.772 

1*066 

1.369 

-2.092 

130.937 

1.124 

1.U87 

-2.448 

90 • U  0  L>' 

0.687 
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TABLE  XIX 


RESULTS  FOR  V  =  U.250.  V/cp  =  4.Q.  PQ  ■  10.0k 


Pt . 

m£ 

1* 

0.000 

1** 

o.ooo 

1*** 

0.000 

2 

-0.263 

3 

-0.527 

4 

-0.791 

5 

-0.263 

6 

—0.527 

7 

-0.791 

8 

-0.527 

9 

-0.791 

10 

-0.791 

-0.527 

12**** 

-0.791 

13 

-0.930 

My 

01/k 

0.000 

-2.593 

0.000 

-9.983 

0.000 

-lu.ooo 

0.000 

-7.682 

0.000 

-6.513 

0.000 

-5.883 

0.050 

-9.339 

0.050 

-7.392 

0.050 

—6.336 

0.100 

-8.738 

0.100 

-7.268 

0.150 

-8.179 

0.076 

-8.149 

0.115 

-7.620 

0.000 

-5.596 

02/k 

O3  /k 

-0.866 

—0 . 866 

-8.251 

-8.251 

-8.234 

-B. 251 

-7.283 

-7.434 

-5.901 

-6.796 

-4.533 

—6  •  296 

-7.607 

-7.607 

—6 .508 

-6.B46 

-5.187 

-6.251 

-7.007 

-7.007 

—5 » 640 

-6.298 

-6*448 

-6.448 

-6.695 

-6.918 

-5.777 

-6.332 

-3.945 

-6.07? 

0° 


75.522 

1.000 

75.522 

1.000 

90.000 

1.014 

90.000 

0.2U1 

89.999 

0.457 

90.000 

0.921 

75.522 

1.000 

127.183 

0.445 

64.304 

0.640 

75.529 

0.999 

133.341 

0.818 

75.528 

0.999 

118.565 

0.783 

128.929 

0.945 

90.000 

1.118 

* 

For 

4P 

s 

4> 

< 

4 

** 

For 

$ 

< 

4> 

< 

4S 

*** 

For 

4s 

< 

4* 

< 

n 

**** 

At 

4s 

At  =  0.043k 
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RESULTS  FOR 


4.0k 


TABLE  XX 

V  =  0.350,  V/cp  •  1.5.  P0  ** 


Pt . 

y? 

py 

01/k 

02/k 

O3  /k 

0° 

n 

1* 

0.000 

0*000 

-3.752 

-2.020 

-2.020 

48.189 

1.000 

1** 

0.000 

0.000 

-4.112 

-2.380 

-2.380 

48.189 

1.000 

0.000 

0.000 

-4.000 

-2.493 

-2.380 

90.000 

0.904 

2 

-0.150 

0.000 

-3.442 

-2.139 

-2.061 

90.000 

U.775 

3 

-0.300 

0.000 

-2.963 

-1.830 

-1.785 

90.000 

0.667 

4 

-0.450 

0.000 

-2.550 

-1.562 

-1.547 

90.000 

0.574 

5 

-0.150 

0.100 

-3.960 

-2«227 

-2.227 

48.189 

1.000 

6 

-0.300 

0.100 

-3»41i 

-1.899 

-1.920 

85.469 

0.866 

7 

-0.450 

0.100 

-2.931 

-1.621 

-1.655 

85.712 

0.746 

8 

-0.300 

0.200 

-3.812 

-2.080 

-2*080 

48.189 

1.000 

9 

-0.450 

0.200 

-3.508 

-1.541 

-1.785 

47.700 

1.072 

10 

-0.450 

0.300 

-3.671 

-1.939 

-1.939 

48.139 

1.000 

1 l*#** 

-0.300 

0.101 

-3.418 

-1.900 

-1.923 

85.407 

0.870 

12**** 

-0.450 

0.152 

-3.170 

-1.634 

-1.720 

83.528 

0.862 

13 

-1.600 

0.000 

-0.803 

-0.349 

-0.511 

90.000 

0.230 

For 

4> 

P 

< 

<p 

< 

$ 

** 

For 

< 

<p 

< 

+s 

*** 

For 

< 

<p 

< 

TT 

#*** 

At 

*+ 

At  «  1.082k 
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TABLE  XXI 


RESULTS  FOR 


0.350.  V/c. 


NOT  TO  SCALE 


Pt. 

y? 

yy 

1* 

0.000 

0.000 

2** 

0.000 

0.000 

1*** 

0.000 

0.000 

2 

-0.150 

0.000 

3 

-0.300 

0.000 

4 

-0.450 

0.000 

5 

-0.150 

0.100 

6 

-0.300 

0.100 

7 

-0.450 

0.100 

8 

-0.300 

0.200 

9 

-0.450 

0.200 

10 

-0.450 

0.300 

11##** 

-0.300 

0.101 

12**»* 

-0.450 

0.152 

13 

-1.500 

0.000 

# 

For 

< 

< 

o 

#* 

For 

< 

<f> 

< 

* 

■ 

*** 

For 

*8 

< 

< 

ir 

****  At  <j>, 


=  138.18 
4>s  =  161.32 
$  *  145.30 


a1/k 

°2/k 

a3/k 

'3.752 

-2.020 

-2.020 

■6.112 

-4.380 

-4.380 

6.000 

-4.493 

-4.380 

5.164 

-3.962 

-3.902 

4.444 

-3.499 

-3*488 

3.825 

-3.096 

-3.130 

5.883 

-4.151 

-4.151 

5.123 

-3.595 

-3.690 

4.403 

-3.178 

-3.293 

5.662 

-3.930 

-3.930 

5.206 

-3.122 

-3.487 

5.450 

-3.718 

-3.718 

5.134 

-3.596 

-3.694 

4.769 

-3.190 

-3.390 

1.376 

-1.335 

-1.657 

48.189 

48.185 

90.000 

90.00U 

90.000 

90.000 

48.189 

83.244 

83.082 

48.189 

47.497 

48.189 

83.165 

80.461 

90.000 


1.000 
1.000 
0.904 
0.711 
0.348 
0.411 
1.000 
0.836 
0.6  lb 
1.000 
1.112 
1.000 
0,861 
0.859 
0.173 


At  =  1.082k 


TABLE  XXII 


RESULTS  FOR  V  =  0.350,  V/c_  =  1.5,  p  =  10.0k 

P  o 


Pt . 

vi  y 

a2/k 

a3/k 

o 

CD 

n 

1* 

0.000 

0.000 

-3.752 

-2.020 

-2.020 

48.189 

l.oou 

1** 

0.000 

0.000 

-10.112 

-8.380 

-8.380 

48.189 

1.000 

1#*# 

0.000 

0.000 

-10.000 

-8.493 

-8.380 

90.000 

0.904 

2 

-0.075 

0.000 

-9.277 

-8.035 

-7.967 

90.000 

0.737 

3 

-0.150 

0. 000 

-8.606 

-7.60/ 

-7.582 

90.000 

0.583 

4 

-0.225 

0.000 

-7.984 

-7.208 

-7.225 

90.000 

0.443 

5 

-0.075 

0.050 

-9.920 

-8.187 

-8.187 

48.189 

l.OJO 

6 

-0.150 

0.050 

-9.234 

-7.713 

-7.782 

84.016 

0.858 

7 

-0.2  .5 

0.050 

-8*565 

-7.305 

-7.405 

83.435 

0.  /OU 

8 

-0.150 

0.100 

-9.730 

-7.997 

-7.99/ 

48.188 

1.000 

9 

-0.225 

0.100 

-9.320 

-7.271 

-7.600 

47.540 

1.100 

10 

-0.225 

0.150 

-9.544 

-7.812 

-7.812 

48.188 

1.000 

1 1***# 

-0.150 

0.050 

-9,243 

-7.7L5 

-7.785 

83.952 

0.863 

12***# 

-0.225 

0.076 

-8.902 

-7.339 

-7.505 

81.296 

U*  658 

13#*** 

-0.300 

0.101 

-8.580 

-6.974 

-7.236 

78.91o 

0.  b5i 

14 

-0.300 

0.000 

-7.407 

-6.830 

-6.893 

90.000 

0.314 

* 

For 

< 

* 

-o-i 

** 

For 

$ 

< 

* 

< 

*s 

*** 

For 

*S 

< 

* 

< 

IT 

#*** 

At 

(J,+ 

At  =  1.032k 
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TABLE  XXIII 


RESULTS  FOR  V  =  0.350,  V/cp  =  2.0,  PQ  =  4.0k 


Pt. 

UC 

uy 

01/k 

02/k 

°3/k 

fl° 

n 

1* 

0.000 

0.000 

-3.752 

-2.020 

-2.020 

60.000 

1.000 

1** 

o.ooo 

0.000 

-4.037 

-2.305 

-2.305 

60.000 

1.000 

1*** 

0.000 

0.000 

-4.000 

-2  •  ?43 

-2.305 

90.000 

0.96  ( 

2 

-0.109 

0.000 

—3.585 

-2.096 

-2.074 

90.000 

0.866 

3 

-0.218 

0.000 

-3.214 

-1.874 

-1.866 

90.000 

o.776 

4 

-0.299 

0.000 

-2.964 

-1.723 

-1.726 

90.000 

0.715 

5 

-0.109 

0.050 

—3 .944 

“2.212 

-2.212 

60.000 

1.000 

6 

-0.218 

0.050 

-3.530 

-1.980 

-1.986 

88.108 

0.893 

7 

-0.299 

0.050 

-3.254 

-1.820 

-1.833 

88.213 

0.823 

8 

-0.218 

0.  100 

-3.852 

-2.120 

-2.120 

59.997 

1.000 

9 

-0.299 

0.100 

-3.638 

-1.861 

-1.954 

58.253 

l.UOO 

10 

-0.299 

0.137 

-3.786 

-2.054 

-2.054 

59.9  n 

1.000 

1 1***# 

-0.218 

0.054 

-3.557 

-1.988 

-1.996 

87.938 

0.903 

12**** 

-0-299 

0.074 

-3.408 

-1.866 

-1.890 

87.21b 

0.883 

13**** 

-1.500 

0.371 

-1.836 

-0.522 

—0 .717 

79*975 

0.708 

14 

-1.500 

0.000 

-0.890 

-0.401 

-0.53  7 

90.000 

u.252 

* 

For 

*P 

< 

4> 

< 

S 

T 

** 

For 

$ 

< 

<f> 

< 

CO 

•e- 

*** 

For 

< 

< 

ir 

**** 

At 

<f>e 

At  =  0.847  k 
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TABLE  XXIV 

RESULTS  FOR  V  =  0.350,  V/cp  =  2.0,  Pq  =  6.0k 


14  4  3  2  i 


Pt„ 

uy 

c^/k 

a2/k 

c3/k 

o 

CD 

n 

1* 

0.000 

0.000 

-3.752 

-2.020 

-2.020 

60. UOO 

1.000 

1** 

0.000 

0.000 

-6.037 

-4.305 

-4.305 

60.000 

1.000 

1*#* 

0.000 

0.000 

-6.000 

-4.343 

—4 .305 

90.000 

0.967 

2 

-0.437 

0.000 

-3.875 

-3.065 

-3.114 

90.000 

L  •  4  5  3 

3 

-0.765 

0.000 

-2.791 

-2.392 

-2.500 

90.000 

0.206 

4 

-1.311 

0.000 

-1.624 

-1.615 

-1.819 

89.999 

0.115 

5 

-0.437 

0.200 

-5.495 

-3.763 

-3.763 

60.000 

1.000 

6 

-0.765 

0.200 

“4.408 

-2.404 

-2*907 

52.029 

1.042 

7 

-1.311 

0*200 

-2.321 

-1.804 

-1 . 966 

81.443 

0.264 

8 

-0.765 

0.350 

-5.120 

-3.388 

-3.388 

59.996 

1.000 

9 

-1.311 

0.350 

-3.592 

-1.470 

-2.182 

49.447 

1.080 

10 

-1.311 

0.600 

-4.550 

-2.818 

-2.818 

59.997 

1.000 

11#**# 

-0.765 

0.189 

-4.023 

-2.685 

-2.879 

79.699 

0.723 

12***# 

-1.311 

0.324 

-3.065 

-1.818 

-2.138 

74.556 

0.647 

13#*#* 

-1.550 

0.383 

-2.728 

-1.508 

-1.869 

72.975 

0.626 

14 

-1.550 

0.000 

-1.383 

-1.269 

-1.614 

90.000 

0.175 

For 

4  <  A 

VP 

<  1 

*# 

For 

-©- 

V 

i-e- 

#**  For  <p„  <  <j>  <  IT 
o 

****  At  4>g 


At  =  o • 847  k 
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TABLE  XXV 


RESULTS  FOR  V  =  0.350,  V/cp  =  2.0,  p  =  10.0k 


9  I 


Pt. 

py 

c^/k 

02/k 

o3/k 

0° 

n 

1* 

o.oou 

o.ouo 

-3.752 

-2.020 

-2.020 

60.UUU 

l.uuo 

1** 

0.000 

0.000 

-10.037 

-8.305 

-8.305 

60.000 

1.000 

1*** 

0.000 

0.000 

-10.000 

-8.343 

-8.305 

90.000 

0.967 

2 

-0.218 

0.000 

-8.036 

-7.171 

-7.208 

90.000 

0.489 

3 

-0.437 

0.000 

-6.458 

-6.213 

-6.320 

90.U0O 

0. 122 

4 

-0.655 

0.000 

-5.430 

-5.190 

-6.602 

90.000 

0 . 207 

5 

-0.218 

0.100 

-9.575 

-7.843 

-7.843 

60.000 

l.ouo 

6 

-0.437 

0.100 

-7.826 

-6.597 

-6.795 

79.099 

0.660 

7 

-0.655 

0.100 

-6.281 

-5.722 

-5.948 

73.468 

0.281 

8 

-0.437 

0.200 

-9.133 

-7.401 

-7.401 

59.994 

1.000 

9 

-0.655 

0.200 

-7.875 

-5.800 

-6.400 

50.897 

1.067 

10 

-0.655 

0.300 

-8.712 

-6.980 

-6.980 

59.994 

1.000 

11*### 

—0 .437 

0.108 

-7.957 

-6.624 

-6.839 

78.643 

0.715 

12###* 

-0.655 

0.162 

-7.139 

-5.867 

-6.216 

73.189 

0.657 

13 

-1.500 

0.000 

-3.489 

-2.225 

-3.885 

90.000 

0. 866 

*  For  <j>p  <  4>  <  ? 

#*  For  ^  <  <p  <  <PS 

**#  For  (j)g  <({><_  it 

....  At 

At  *  0.847k 


67 


TABLE  XXVI 


RESULTS  FOR  V  =  0.350,  V/cp  =  4.0,  pQ  =  ^»0k 


Pt. 

1* 

0.000 

1** 

0.000 

0.000 

2 

-0.470 

3 

-0.940 

4 

-1.410 

5 

-0.470 

6 

-0.940 

7 

-1.410 

8 

-0.940 

9 

-1.410 

10 

-1.410 

-0.940 

12**** 

-1.410 

yy 


0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.100 

0.100 

0.100 

0.200 

0.200 

0.300 

0.113 

0.170 


°l/k 


-3. 752 
-5.005 
-5.U0U 
-3.124 
-1.952 
-1.218 
-4.583 
-2.840 
-1.804 
-4,194 
-2.534 
-3.838 
-2.993 
-2.339 


°2/k 


-2.020 

-3.273 

-3.278 

-2.189 

-1.483 

-1.016 

-2.851 

-1.869 

-1.180 

-2.463 

-1.620 

-2.106 

-1.938 

-1.418 


C3/k 


-2.020 

-3.273 

-3.273 

-2.235 

-1.578 

-1.158 

-2.851 

-1.897 

-1.294 

-2.463 

-1.587 

-2.106 

-1.958 

-1.481 


75.522 

75.522 

90.000 

90.000 

90.000 

90.000 

75.522 

98.950 

110.612 

75.521 

73.550 

75.521 

98.376 

106.188 


# 

For 

*P 

< 

* 

< 

<J> 

#* 

For 

$ 

< 

4> 

< 

*s 

*** 

For 

*S 

< 

<t> 

< 

ir 

**** 

At 

♦; 

At 

*  0 

•  431k 

n 


1.000 

1.000 

0.995 

0.526 

0.247 

0.103 

1.000 

0.553 

0.332 

0.999 

0.537 

0.999 

0.603 

0.514 
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TABLE  XXVII 


pt.  u?  py  a]/k  °2^k 


1*  0*000  0*000  -3*752  -2*020 

1**  0.000  0.000  -7.005  -5.273 

1*#*  0.000  0.000  -7.000  -5.278 

2  “0.470  0.000  -4.374  -3.754 

3  -0.940  0.000  -2.765  -2.732 

4  -1.410  0.000  -2.112  -1.705 

5  -0.470  0.100  -6.414  -4.682 

6  -0.940  0.100  -4.009  -3.273 

7  -1.410  0.100  -2.722  -2.145 

8  -0.940  0.200  -5.870  -4.139 

9  -1.410  0.200  -3.547  -2.958 

10  -1.410  0.300  -5.371  -3.639 

11#***  -0.940  0tll3  -4.217  -3.376 

12****  -1.410  0.170  -3.373  -2.575 


°3/k 


-2.020 

-5.273 

-5.273 

-3.C2U 

-2.900 

-2.311 

-4.682 

-3.347 

-2.502 

-4.139 

-2.912 

-3.639 

-3.433 

-2*764 


G° 


75.522 

75.522 

90.000 

90.000 

89.999 

89.999 

75.522 

107.300 

46.569 

75.521 

71.222 

75.521 

105.202 

119.950 


n 


1.000 

1.000 

0.995 

0.340 

0.088 

0.309 

1.000 

0.405 

0.290 

0.999 

0.354 

0.999 

0.470 

0.417 


#  For  <j)p  '  4>  <  $ 

**  For  $  <  4>  <  4>g 

##*  For  4>s  < 

*#**  At  4>g 

At  *  0.431k 
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TABLE  XXVIII 


RESULTS  FOR  V  =  0.350,  V/cp  =  4.0,  pQ  =  10.0k 


4  3  2  1 


rP 


Pt . 

US 

uy 

°l/k 

1* 

0.000 

0.000 

-3.752 

i** 

0.000 

0.000 

-10.005 

1*** 

0.000 

O.OUO 

-10,000 

2 

-0.470 

0.000 

-6.248 

3 

-0.940 

0.000 

-4,688 

4 

-1.500 

0.000 

-3.615 

5 

-0.470 

0.  loo 

-9.161 

6 

-0*940 

0.100 

-5.876 

7 

-1.500 

0.100 

-4.219 

8 

-0.940 

0.200 

-8.384 

9 

-1.500 

0.200 

-4.761 

10 

-1.500 

0.319 

-7.542 

1 1**** 

-0.940 

0.113 

-6.131 

12**** 

-1.500 

0.181 

-4.968 

•»r 

For 

< 

4> 

< 

$ 

** 

For 

-e-i 

< 

4> 

< 

s 

*#* 

For 

+s 

< 

<P 

< 

ir 

*#** 

At 

♦s 

°2/k 

°3/k 

o 

CD 

a 

-2.020 

-2.020 

75.522 

l.uoo 

-8.273 

-8.273 

75.522 

1.000 

-8.278 

-8.273 

90.000 

0.995 

-6.101 

-6.198 

90.000 

0.0  75 

-3.904 

’4.083 

90.000 

0.5,18 

-2.225 

3.920 

90.000 

0.903 

-7.429 

7.429 

75.522 

1.000 

-5.265 

-5021 

128.944 

0 . 306 

-2.975 

-4.140 

68.375 

0,690 

-6,652 

-6.652 

75  *  520 

0.999 

-4.552 

—4.649 

86.927 

u.  104 

-5.811 

-5.811 

75.521 

0.999 

-5.454 

-5 . 644 

121.934 

0.349 

-3.880 

-4.528 

51.370 

0.547 

Ax  =  0.431k 


case  the  value  0  stays  at  60.00°  for  <f>  <  c|)  <  and  at 

Jt  j 

90.00°  for  $  <  <j>  <  ti. 

Looking  at  results  for  other  grid  points  it  is  seen 
that  at  point  10  in  Table  XIV  the  plastic  shock  has  reached 
a  small  tensile  value,  indicating  that  the  requirement  of 
compressive  shocks  is  violated  from  there  on. 

For  the  higher  loads  PQ/k  =  6.0  and  10.0,  Tables  XV,  XVI, 
the  shock  remains  compressive  in  the  entire  range  of  the  tables. 
Table  XV  for  PQ/k  =  6.0  is  limited  only  by  the  condition 
p£  -1.7  due  to  the  truncation  of  the  expansions.  For  the 
highest  load  PQ/k  =  10.0,  the  validity  is  limited  to  p£  -1.0 
because  the  yield  condition  on  the  surface,  beyond  point  14 
would  be  violated  excessively. 

For  other  values  of  V  and  V/c_,  the  situation  is  generally 

sr 

similar.  Violation  of  the  yield  condition  for  low  values  of 

PQ/k  occurs  first  on  the  S-front,  while  for  higher  values  of 

p  /k  the  violation  appears  first  on  the  surface, 
o 

2 .  Range  of  depth  for  which  results  apply 

The  previous  discussions  considered  the  validity  of  results 
with  respect  to  a  nondimensional  depth  py.  It  is  appropriate 
to  consider  values  of  p  of  physical  interest  to  find  out  for 
what  actual  depth  y  results  have  been  obtained. 
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800  ft. 


Let  t  =  0.2  sec,  V  =  4000  ft/sec,  then  -jj  =  VtQ  = 

For  the  typical  values  V  =  0.25  and  V/cp  =  2.0  discussed  above, 

information  on  peak  pressures  for  values  of  the  surface  load 

p  =  4.0k  and  6.0k  are  contained  in  Tables  XIV  and  XV  up  to 
o 

about  py  =  0.6  or  for  a  depth  of  480  feet.  A  time  history 

can  be  plotted  using  interpolation  for  the  depths  p  y  =  0.5 

or  400  feet.  For  the  same  values  of  V  and  V/cp  ,  but 

p  /k  =  10.0,  the  limit  of  applicability  is  py  =  0.3,  or  about 
ro 

240  feet. 

In  general,  for  larger  values  of  V/cp  and  for  larger 

values  of  p  /k  the  depth  for  which  the  results  apply  decreases, 
o 


SECTION  V 


CONCLUSIONS 

An  analytical  approach  has  been  presented  to  describe 
the  effect  of  a  decaying  pressure  moving  with  superseismic 
velocity  V  _>  1.5  cp  on  the  surface  of  a  half-space  of  an 
elastic-plastic  material  subject  to  the  von  Mises  yield  con¬ 
dition.  Numerical  results  have  been  obtained  for  a  wide 
range  of  surface  pressure  and  material  properties.  Due  to 
simplifications  in  the  analysis  the  results  obtained  are 
approximate  and  restricted  to  target  points  of  limited  depth, 
as  discussed  in  Section  IV. 

Information  on  the  stress  field  for  exponential  decay 
of  the  applied  load  is  contained  in  Tables  II  -  XXVIII  and 
in  the  accompanying  figures.  Comparing  the  effects  of  step 
loads  given  in  Ref.  [1]  with  the  new  results  for  decaying 
loads,  it  is  found  that  in  the  latter  case  the  peak  and  sub¬ 
sequent  pressures  decrease  with  depth,  which  is  not  the  case 
for  a  step  load.  However,  the  decrease  found  at  any  depth  y 
is  smaller  than  the  decrease  of  the  applied  pressure  at  the 
point  on  the  surface  directly  above  the  target  point,  i.e., 
at  a  distance  equal  to  (-y/tg  <{>)  behind  the  shock  front.  This 
is  graphically  shown  in  Figs.  10  to  18.  While  the  attenuation 
is  therefore  less  important  than  one  might  expect,  it  is  still 
greater  than  found  in  a  similar  analysis,  Ref.  [3],  for  an 
elastic-plastic  material  of  the  Coulomb  type. 
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APPENDIX  I  -  Special  Case,  m  =  0 


Although  an  approach  analogous  to  that  used  in  the  text 
for  m  >  0  may  be  used  for  m  =  0,  it  is  considerably  simpler 
to  derive  expressions  for  a  step  pressure  using  the  simple 
configuration  shown  in  Fig.  8,  consisting  of  three  shock 
fronts  only. 


For 


v  Ao 
l-\>  k 


where 


/T(l-v) 

1-  2  V 


For  <{>  <  <j>  <  4>( 


AaP  ,  Ao 

T"  +  T 


\a  ,  v_ .  ,  Aa 
IT  +  T 


where 


,  /3  . 2(l-v)  sin  2<t>?  sin  2 

2  C  1- 2v  t0S  2^S  +  cos  26  _ 


cos  2<J>, 
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I  ” 


APPENDIX  II  -  Elastic  Solutions 


The  expressions  for  the  stresses  in  a  purely  elastic 


material  subjected  to  a  step  pressure  may  be  found  in  Ref.  [1] 
The  equivalent  expressions  for  a  decaying  surface  pressure  are 


For  <j>p  <  <P  <  <PS 


k  =  T  cos  24>s  e 


p(5-?d) 


°x  A  o'  f  ,  .  2,  2.  ,  P(5"5P* 

—  =  (1-2  sin  s  cot  4>p]  e 


°z  AoP  rl  ,l-2v. 

T  =  IT  [2  sin  2+Pl  e 


T  A0P  r„  2.  .  . 

-  -jj-  [2  sin  <f>s  cot  <j>pJ  e 


Pq  (1-v)  ,  cos  2<J) g 

k  [cos2<j>s  +  ( l-2v)  cos2(<j>s  -  4>p) 


(1-v)  ] 


y  cos  <J>. 

£  =  - f 

P  sin 


(100) 


For  4>s  <  (j)  <  tt 


o  p  _ 
_Z  =  _  12  ey? 
k  k 


(101) 
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(102) 


?Z  ■  Aap  ( l~2v)  4  ^ 

k  k  2  ( l-v)  sin  2<^P  e  (103) 


f  2  sin2<f>g  cot  4>p]e 


M(C-5p) 


cos  2<j>s  e 


m(€-Cs) 


AtS  =  1- 2v  sin  2(^p  AoP 
k  “  2  (1-v)  cos  2<J>  k 

u 

y  cos  <j>c 

r  =  _ s 

o  sin  tf> 

O 


(104) 

(105) 

(106) 
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